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^ The trace anomaly of quantum fields in electromagnetic or gravitational backgrounds im- 

I tensor. Considering first the axial anomaly and using QED as an example, we compute 

the full one-loop triangle amplitude of the fermionic stress tensor with two current vertices, 
(T^'^J^J'^), and exhibit the scalar pole in this amplitude associated with the trace anomaly, 
in the limit of zero electron mass m 0. To emphasize the infrared aspect of the anomaly, 
we use a dispersive approach and show that this amplitude and the existence of the mass- 
less scalar pole is determined completely by its ultraviolet finite terms, together with the 
^ requirements of Poincare invariance of the vacuum, Bose symmetry under interchange of 

J" and J^, and vector current and stress tensor conservation. We derive a sum rule for 
the appropriate positive spectral function corresponding to the discontinuity of the triangle 
amplitude, showing that it becomes proportional to S{k^) and therefore contains a mass- 
less scalar intermediate state in the conformal limit of zero electron mass. The effective 
action corresponding to the trace of the triangle amplitude can be expressed in local form by 
the introduction of two scalar auxiliary fields which satisfy massless wave equations. These 
massless scalar degrees of freedom couple to classical sources, contribute to gravitational 
scattering processes, and can have long range gravitational effects. 
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I. INTRODUCTION 

Quantum effects are most commonly associated witli sliort distance pliysics. Tlie basic reason for this 
is that the fluctuations of most fields have a finite correlation length, and hence their correlations fall off 
rapidly at large distances. For relativistic fields which are massive, the correlation length is m"-'^, and the 
fall off is exponential, e~'^^. This is an elementary example of decoupling. In the limit of very large mass, 
m — > oo, the quantum effects of a heavy field become negligible at any finite distance scale. 

In the opposite limit of massless fields m ^ 0, the correlation length becomes infinite. Decoupling no 
longer holds, and it becomes possible for quantum correlations to extend over very great distances on even 
macroscopic scales. Such infrared effects are more pronounced the lower the spacetime dimensionality. 
In d = 2 the two-point propagator function of a free massless scalar field, rather than falling off, grows 
logarithmically in the separation of the points. Massless conformal field theories in d = 2 have been 
studied extensively by a variety of methods, and it is clear that their fluctuations are non-negligible and 
have important physical consequences in the infrared [1]. 

In massless field theories in two dimensions, the conformal group algebra and its central extension play 
important roles. The central term, a Schwinger term in the commutation algebra of stress-tensors in flat 
space, may also be recognized as the trace anomaly of the stress-energy tensor defined in curved space 
[U [2] . The anomaly in the trace of the stress-energy tensor, corresponds to a well-defined additional term 
in the effective action which has long range effects [3 . As one illustrative example of these infrared effects, 
one can show that the curved space anomaly modifies the critical scaling exponents of the two-dimensional 
Ising model at its second order phase transition point This modification of the critical exponents is 
associated with the fluctuations of the spacetime metric at large distance scales. The gravitational metric 
fluctuations may be described by an additional massless scalar field, the conformal or Liouville mode, 
whose dynamics is generated and required by the conformal anomaly. 

In two dimensions the central term related to the anomaly in curved space can be seen already in 
the two-point correlation of stress tensors, (T^''(a;)T"^(y)). The meaning of this term is most clearly 
deduced from the momentum space representation, where the corresponding amplitude exhibits a massless 
pole in the conformal limit [5 . This massless pole corresponds to a 6{k'^) in the corresponding imaginary 
part, describing a propagating massless scalar degree of freedom in the two-particle intermediate state of 
the cut diagram. The same kinematics applies to the pole and massless scalar state in the correlator of 
electromagnetic currents {J^{x)J'^{y)) in the Schwinger model of two-dimensional massless electrodynamics 
[HI E]. The infrared effects of the anomaly may be understood as the result of the fluctuations of this 
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additional massless degree of freedom. 

In dimensions greater than two, infrared effects due to anomalies are both more subtle, and somewhat 
less well studied. Since the correlation functions of canonical free theories fall off as power laws at large 
distances for d > 2, at first sight there would seem to be little possibility of enhanced infrared effects in 
higher dimensions. QCD is a notable counterexample, where the growth of the effective coupling at larger 
distances leads to large quantum fluctuations and infrared confinement. The renormalization group flow 
f3{g'^) of the coupling arises from the same breaking of scale invariance by quantum fluctuations which give 
rise to the conformal anomaly pi . 

In the most familiar case of the axial anomaly, a massless pseudoscalar pole does appear in the triangle 
amplitude (J^ J"J^) [9l |T0], in the chiral limit of vanishing fermion mass, a feature we review in the 
next section. This example of the axial anomaly in massless quantum electrodynamics (QED) shows that 
infrared relevant fluctuations due to anomalies can occur in d = 4, and that triangle amplitudes are the 
simplest ones to reveal these effects. In QCD the lightest pseudoscalar state is the pion, whose mass vanishes 
in the chiral limit of zero quark mass. By identifying this state with the massless pole appearing in the 
perturbative anomaly in the chiral limit, the low energy rate of neutral pion decay, vr*^ 27 is determined 
by the short distance colored quark degrees of freedom in the one- loop {J^J'^P) amplitude [mill]. The 
agreement of the measured rate with the coefficient obtained with Nc = 3 quarks is a striking confirmation 
of both QCD and the infrared effects of the anomaly. This well-known example of anomaly matching |13j 
shows that anomalies can provide a mechanism for short distance quantum degrees of freedom to have long 
distance or low energy consequences. 

Although the special role of the triangle diagram in d = 4 has been emphasized in [U] in the context 
of the chiral anomaly some time ago, to date there has been no clear indication of a massless pole or 
infrared degrees of freedom in flat space amplitudes involving the energy-momentum tensor. It is known 
that in d = 4 the trace anomaly in curved space involves geometric invariants that are quadratic in the 
Riemnann curvature tensor [2, 15j. This has the immediate consequence that the simplest amplitude in 
four dimensional flat spacetime that can show any direct evidence of the full curved space anomaly is the 
three-point function of stress tensors, {T^^'^ {x)T°'^ {y)T'^^ (z)) , indicating again the importance of triangle 
amplitudes in d = 4. 

In this paper we address the possibility for low energy quantum effects in gravity, analogous to those 
in gauge theories, due to the corresponding trace anomaly, and in particular for additional massless scalar 
degrees of freedom with long range effects which can modify the predictions of classical General Relativity 
on macroscopic and even cosmological scales. We present a complete calculation of the one-loop triangle 
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amplitude {T^'^ J'^Jl^) in QED, for all values of the kinematical invariants. This amplitude contains the 
same basic kinematics as both the more familiar chiral triangle (Jg J"J^), and the more complicated 
amplitude (TTT) involving three stress tensors. The {TJJ) amplitude is sensitive to the trace anomaly 
of the one-loop stress tensor expectation value in a background electromagnetic potential (rather than 
a gravitational background curvature). By calculating this amplitude for arbitrary electron mass m, both 
the decoupling limit m ^ oo, and the conformal limit m — > 0, where the massless pole in the amplitude 
appears can be studied. It is the latter limit that reveals the consequences for low energy gravity. 

Following methods that have been used previously for the chiral anomaly |16| [T7] , we show that the 
(TJJ) amplitude can be determined completely from general principles of Poincare invariance of the 
vacuum, Bose symmetry under interchange of J° and J^, and the Ward identities of vector current and 
stress-energy conservation, once its finite tensor components are given. These finite components can be 
determined unambiguously from the imaginary part of the cut triangle amplitude, with the real part 
obtained by dispersion relations which require no subtractions (other than charge renormalization in one 
particular component). This dispersive approach based upon the finite parts of the amplitude emphasizes 
the infrared aspect of the anomaly, making it clear that the anomaly is finite, well defined and uniquely 
determined, independent of UV regularization scheme, provided only that the amplitude is defined in a 
way consistent with the non-anomalous low energy symmetries of the theory. 

In the conformal limit of massless QED, the two-particle intermediate state of the cut triangle diagram 
has a delta function contribution at k'^ = 0. Because this state couples to the stress tensor, it contributes 
to gravitational scattering amplitudes at arbitrarily low energies. We demonstrate that the trace part of 
the (TJJ) amplitude containing this massless intermediate scalar state and its gravitational couplings may 
be described by the introduction of local massless scalar degrees of freedom, which render the trace part 
of the one-loop effective action local. The auxiliary field description introduced recently in the context of 
curved space |I8] reproduces the trace part of the amplitude exactly in flat space, and the massless pole in 
the trace part of the flat space (TJJ) amplitude is precisely the propagator of these scalar fields. 

In QED the scalar state may also be understood as a two-particle correlation of e^e~ which in the 
massless limit move coUinearly at the speed of light in a total spin-0 configuration. When the electron 
mass is non-zero, the singularity at A;^ = is replaced by a resonance with a width of order m^. However, 
the corresponding spectral function obeys a sum rule, which shows that although broadened, and eventually 
decoupled for larger m, the scalar state survives deformations away from the conformal limit. In this sense 
it behaves analogously to the pion in QCD. 

The paper is organized as follows. In the next section we review the axial anomaly in QED in four 
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dimensions, using the spectral representation and dispersion relations to emphasize its infrared character, 
exhibiting the massless 0~ intermediate state, and the finite spectral sum rule in this case. In section 3 we 
give the auxiliary field description of the chiral amplitude, showing that the massless pseudoscalar state 
can be described by a local effective field theory. In section 4 we turn to the main task of evaluating the 
(TJJ) amplitude in QED. Imposing the Ward identities, we show that the full amplitude is determined 
by its finite terms and imaginary parts for any m and its three kinematic invariants, independently of 
any specific UV regularization method. In section 5 we evaluate its trace, isolate the anomaly and discuss 
its relation to the /3 function and scaling violation. In section 6 we give the spectral representation of 
the {TJJ) amplitude, derive the corresponding finite sum rule, and show that a 6{k'^) appears in the 
appropriate spectral function in the conformal limit of massless electrons. In section 7 the foregoing results 
are compared with the auxiliary field representation of the anomaly given in [18, , and shown to coincide 
exactly in the trace sector. In section 8 we show how the anomalous amplitude contributes to gravitational 
scattering of photons by a source, prove that the anomaly pole induces a massless scalar interaction and 
propagating intermediate state in this scattering process, and provide the effective action description of 
the scattering by scalar exchange. Finally, section 9 contains a concise summary of our results. Technical 
details of extracting the finite parts of the {TJJ) amplitude are given in Appendix A, while the proofs of 
some identities needed in the text are given in Appendix B. 



II. THE AXIAL ANOMALY IN QED 

In order to exhibit the relationship between anomalies and massless degrees of freedom, we review first 
the familiar case of the axial anomaly in QED in this section |1H \T9\ [20] . Although the triangle anomaly 
has been known for quite some time, the general behavior of the amplitude off the photon mass shell, its 
spectral representation, the appearance of a massless pseudoscalar pole, and its infrared aspects generally 
have received only limited attention |9l |2T] . It is this generally less emphasized infrared character of the 
axial anomaly upon which we focus here. 

The vector and axial currents in QED are defined by^ 

J^'{x) = V;(x)7^V(^) • (2.1a) 
J^{x) = i){x)-f^'-^^^{x) . (2.1b) 

^ We use the conventions that {7'', 7"} = —2 3^" = 2diag (H ), so that 7" = (7°)^ and 7^ = i7''7^7^7'' = (7^)''" are 

hermitian, and tr(7^7'^7'^7''7'^) = — 4ie^'''"', where e'"''"' = —e^^pa is the fully anti-symmetric Levi-Civita tensor, with 
£0123 ~ +1- 
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The Dirac eq., 

- i7^(5/, - ieA^)7p + m^j = . (2.2) 
implies that the vector current is conserved, 

9^J'^ = 0, (2.3) 

while the axial current apparently obeys 

df,J^ = 2imi}-f^ij (classically). (2.4) 

In the limit of vanishing fermion mass m — > 0, the classical Lagrangian has a f/c/i(l) global symmetry 
under ip — > e^°''^^ip, in addition to U{1) local gauge invariance, and Jg is the Noether current corresponding 
to this chiral symmetry. As is well known, both symmetries cannot be maintained simultaneously at the 
quantum level. Let us denote by {Jt^{z))^ the expectation value of the chiral current in the presence of 
a background electromagnetic potential A^. Enforcing f7(l) gauge invariance (2.3) on the full quantum 
theory leads necessarily to a finite axial current anomaly, 

in a background electromagnetic field. 

The second variation of {Jt^(z))^ with Afj_ then set to zero. 



r^-/^(.:x,,)^-. 

6Aa{x)dAf3{y) 



-^{ie)\TJ^{z)^[x)JP{y))\.^ (2.6) 



A=0 

A=0 



is thus the primary quantity of interest. By translational invariance of the Minkowski vacuum at ^ = 0, this 
amplitude depends only upon the coordinate differences x — z and y — z. Hence with no loss of generality 



we may fix z = 0. Taking the Fourier transform of (2.6) and removing the factor of total momentum 
conservation, (27r)'^(5^(A; — p — q), we obtain 



Jp-x+iq-y ^ {•^5 (Q))/ 



d^x / d^'ye^^-^+^i-y {T4{Q)r{x)J^{y))\.. . (2.7) 



5A^{x)5Ap{y) 



At the lowest one-loop order it is given by the triangle diagram of Fig. [T] plus the Bose symmetrized 
diagram with the photon legs interchanged. The chiral current expectation value in position space can be 
reconstructed from this momentum space amplitude by 
i f d'^p f d^q 



(-^rC^)). = ^ y y / '^'^ / A e-^^-(^-^) e-^-(^-^) T^'^^p, q) A^{x)A,{y) + ... {21 
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FIG. 1: The Axial Anomaly Triangle Diagram 

up to second order in the gauge field background Afj,. 

Elementary power counting indicates that the triangle diagram of Fig. [T] is superficially linearly di- 



vergent. The formal reason why (2.3) and (2.4) cannot both be maintained at the quantum level is that 



verifying them requires the ability to shift the loop momentum integration variable / in the triangle am- 
plitude. Because the diagram is linearly divergent, such shifts are inherently ambiguous, and can generate 
finite extra terms. It turns out that there is no choice for removing the ambiguity which satisfies both the 
vector and chiral Ward identities simultaneously, and one is forced to choose between them. Thus although 
the ambiguity results in a well-defined finite term, the axial anomaly has most often been presented as in- 
herently a problem of regularization of an apparently ultraviolet linearly divergent loop integral |1H [T9\ [20] . 

There is an alternative derivation of the axial anomaly that emphasizes instead its infrared character. 
The idea of this approach is to use the tensor structure of the triangle amplitude to extract its well-defined 
ultraviolet finite parts, which are homogeneous of degree three in the external momenta p and q. Then the 
remaining parts of the full amplitude may be determined by the joint requirements of Lorentz covariance, 
Bose symmetry under interchange of the two photon legs, and electromagnetic current conservation. 



(2.9) 



at the two vector vertices. By this method the full one-loop triangle contribution to r^"^(p, g), becomes 
completely determined in terms of well-defined ultraviolet finite integrals which require no further regu- 
larization |16l [T7] . The divergence of the axial current may then be computed unambiguously, and one 



obtains (2.5) in the limit of vanishing fermion mass [TT]. There is of course no contradiction between these 
two points of view, since it is the same Ward identities which are imposed in either method, and in the 
conformal limit of vanishing fermion mass the infrared and ultraviolet behavior of the triangle amplitude 
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TABLE I: The 8 three index psuedo-tensor monomials into which r^"''(p,g) can be expanded 
are one and the same. 

Since we will apply this method to the trace anomaly amplitude in the next section, let us first review 
the calculation in the axial current case. One first uses the Poincare invariance of the vacuum to assert that 
Tf^°'^{p, q) can be expanded in the set of all three-index pseudotensors constructible from the p and q, with 
the correct Lorentz transformation properties. There are exactly eight such pseudotensors, which are listed 
in Table |lj the first two of which are linear in p or q, while the remaining six are homogeneous of degree 
three in the external momenta. Since the amplitude ^^"'^{p, q) has mass dimension one, any regularization 
ambiguity can appear only in the coefficients of the tensors which are linear in momenta, i.e. e'^°'^^p\ 
and The coefficients of these tensors have mass dimension zero and are therefore potentially 

logarithmically divergent. On the other hand, the remaining six tensors in Table |lj homogeneous of degree 
three in p and q can appear in T^^°^^{p,q) multiplied only by scalar loop integrals with negative mass 
dimension, —2, which are completely convergent in the ultraviolet. If these scalar coefficient functions can 
be extracted unambiguously, then vector current conservation can be used to determine the coefficients of 



the remaining two tensors of dimension one. Indeed the general amplitude satisfying (2.9) must be a linear 



combination of only the six linear combinations defined below and listed in Table llTl Since the tensors 



^fiafSXp^ and e^'^^^qx appear only in those linear combinations which satisfy (2.9), their coefficients are 



determined unambiguously by the finite coefficients multiplying the tensors of degree three. 



i 




1 


~p ■ q e^^'^Ap^ _ pi3 


2 




3 


p'' v°'^{p, q) 


4 


p.qet'°'P^qx + q°' v^''^{p, q) 


5 


_q2^pal3Xp^ _ qPyt^a{p^q) 


6 


qt^ 



TABLE IL The 6 third rank pseudotensors obeying (2.121 
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To make the procedure of extraction of finite terms of the amplitude completely unambiguous, one may 
first calculate the imaginary part of the cut triangle amplitude in Fig. [2] at timelike fc^, which is finite, and 
then construct the real part by a dispersion relation. For the mass dimension —2 terms, the dispersion 
relations constructing the real parts of the amplitude from its imaginary parts are finite and require no 
subtractions O |21] . 



a P 
/WWW 



•^5 




FIG. 2: Discontinuity or Imaginary Part of the Triangle Diagram obtained by cutting two lines 



To construct the tensors satisfying (2.9), let us define first the two index pseudotensor. 



which satisfies 



PaV°''^{p,q) = = q(3v'^'^{p,q) 



(2.10) 

(2.11a) 
(2.11b) 



By taking general linear combinations of the eight pseudotensors in Table |Tj we find then that there are 
exactly six third rank pseudotensors, T^'^^ip, g), z = 1, . . . , 6 which can be constructed from them to satisfy 



the conditions (2.9), 



PaT^'^'^ip, q) = = Tf'°"^{p, q) 



0. 



i = l,...,6, 



(2.12) 



given in Table |II[ Hence we may express the amplitude (2.7) satisfying (2.9) in the form, 

6 



T^^'-P{p,q) = ^f,rr''{p,q), 



(2.13) 



i=l 



where fi = fi{k'^;p'^, q^) are dimension —2 scalar functions of the three invariants, p^ , q^, and /c'^. 
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We note also that the full amplitude (2.7) must be Bose symmetric, 



(2.14) 



Since T-T^{p,q) = T^^'^{q,p) for i = 1,2,3, it follows that the six scalar coefficient functions fi also fall 



into three Bose conjugate pairs, i.e. 

f2[k';p\q^) = Uk';q\p'), 
Hk';p\q^) = Uk';q\p'), 

related by interchange of and q^. 

Actually, owing to the algebraic identity obeyed by the e symbol. 



gOtP^Hvpa _|_ gafi^up(Tl3 _|_ gau^paf3p _|_ gap^a/Sfiu _|_ gOa^/Spup _ q 



(2.15a) 
(2.15b) 
(2.15c) 



(2.16) 



in four dimensions, the six tensors r, are not linearly independent, and form an overcomplete basis. The 



identity (2.16) leads to the relations. 



{p, q) = T^'^iP, q) + T^"'\p, 



paf3f 



pais 



iP,Q) 



pafS 



{p, q) + Tr^ijp,q) 



(2.17a) 
(2.17b) 



Thus the tensors T3 and rg could be eliminated completely by means of (2.17), and the full amplitude 



expressed entirely in terms of the linearly independent and complete basis set of only the four tensors, 
T"2, '^4, and T5. Indeed this has been the general practice in the literature on the axial anomaly |lll [T6| [T9] 
Eliminating these or any other two tensors is not necessary for our purposes, and we choose instead to 
work with the overcomplete set of six tensors listed in Table |ITj This will have the consequence that the 
coefficient functions fi are determined only up to the freedom to choose arbitrary coefficients of the linear 



combinations (2.17), i.e. to shift each of the coefficients fi by an arbitrary scalar function h via the rule. 



fl{k^;p\q')^fl{k';p\q^) + h{k^;p\q') 
f2{k^;p^q^) ^ f2{k^;p^q^) + Kk^;p^q^) 
f3ik^;p^q^) ^ f3{k^;p\q'') - Kk'-p^q^) 



(2.18a) 
(2.18b) 
(2.18c) 



with the shift in fi, f^, fe determined by (2.15) by interchange of p and q . The arbitrary function h drops 



out of the final amplitude by use of (2.17) 
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The computation of the finite coefficients given in the hterature |1H [T6] amounts to a specific choice of 
the arbitrary function h (by the order of the 7 matrices when the trace is performed) , and yields 



/l = /4 



/2 



/a = h 



^ Jo 



dx 



dx 







dx 

l~x 



l-x 



dy 



dy 



dy 



xy 
D 

x(l — x) 
D 

y(i - y) 



(2.19a) 
(2.19b) 



D 



(2.19c) 

0, (2.19d) 
where the denominator of the Feynman parameter integral is given by 

D = p^x(l — x) + q^y{l — y) +2p ■ qxy + rr? = [j? x + y){l — x — y) + xy k"^ + rr? , (2.20) 
strictly positive for rr? > 0, and spacelike momenta, k'^,p'^,q'^ > 0. Thus each of the dimension —2 scalar 



coefficient functions fi in (2.19) are finite, and free of any UV regularization ambiguities, and the full 
amplitude T^'^^{p,q) satisfying 

(i) Lorentz invariance of the vacuum , 



(ii) Bose symmetry (2.14) 



(iii) vector current conservation (2.9) , 

(iv) unsubtracted dispersion relation of real and imaginary parts , 

with the finite imaginary parts determined by the cut triangle diagram of Fig. [2| is given by (2.13) and 



(2.19), without any need of regularization of ultraviolet divergent loop integrals at any step. 

Contraction of the finite amplitude T^'^^{p, q) with the momentum = (p + g)^ entering at the axial 
vector vertex can now be computed unambiguously, and we obtain 



with 



k^T^'''P{p,q)=Av''f{p,q), 



A = p ■ q fi + p^f2 + {p^ + P ■ q)h + p ■ q f4 + q^ h + {p-q + q^)k 
= 2p-qh+p'f2 + q^h, 



(2.21) 



(2.22) 



by (2.19). Substituting the explicit Feynman parameter integrals of (2.19) in (2.22), and using (2.20), 



(2.22 ) becomes 



A{k^;p\q') 

- ^_ 
~ 27r2 ■ 



dx 
1 



, D-m^ 
dy 







dx 







dy 



D 
1 

D 



(2.23) 
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The second term proportional to is what would be expected from the naive axial vector divergence 



(2.4) |T7|. The first term in (2.23) in which the denominator D is cancelled in the numerator is 



dx 



l-x 



dy 



e 

2^' 



(2.24) 



and which remains finite and non-zero in the limit m — > is the axial anomaly. 

Thus the finite anomalous term is unambiguously determined by our four requirements above, and 
may be clearly identified even for finite m, when the chiral symmetry is broken. This construction of the 
amplitude from only symmetry principles and its finite parts may be regarded as a proof that the same 
finite axial anomaly must arise in any regularization of the original triangle amplitude which respects these 
symmetries and leaves the finite parts unchanged. Explicit calculations in dimensional regularization and 
Pauli-Villars regularization schemes, which respect these symmetries confirm this ^22j . 

The spectral representations for the triangle amplitude functions, 

P^(s;P^9^) 



ds 



(2.25) 



A;2 + s ' 

used to compute the finite parts also aid in the physical interpretation of the infrared aspect of the anomaly. 
If one defines the function, 

(p^x + g^y)(l — X — y) -\- m? D 



S{x,y;p^,q^) = 
and substitutes the identity, 

1 

D 



xy 



ds 



xy 



xy{k'^ + s) 



5{s - S) 



(2.26) 



(2.27) 



valid for p ,q ,m > 0, into the expressions (2.19), interchanging the order of the s and x^y integrations. 



the spectral representation (2.25) of the amplitude is obtained, with 

"1 ri-x 



Pi{s;p ,q ) 







P2{s;p ,q ) = — 



P5{s;p ,q ) 



dx 



dx 



dx 



dy 5{s — S) 







l-x 



dy 



1 — X 



l-x 



dy 



1-2/ 



5{s - S) 



5{s - S)= P2{s;q ,p ) 



From the definition of the function S{x, y;p , q ) in (2.26) it follows that 



r, 1 — X r, 1 — y D — m 
2p-q+p^ + q^ ^ = 

y X xy 



k^ + S 



m 
xy 



and therefore, from (2.28), 



2p- qpi+p^ P2 + q^ P5 = {k^ + s)p^ - m^po 



(2.28a) 
(2.28b) 
(2.28c) 

(2.29) 

(2.30) 



13 



where 



PaV^'^p ,q ) = Pi{s;p ,q ) 



■7r2 , 



dx 



1-x 



and 

dy — d{s — S) . 

xy 



(2.31a) 
(2.31b) 



The relation (2.30) with (2.31) for the imaginary part of the cut triangle amplitude can be compared to 
(2.22 )-(2.23) for the corresponding real part. Defined for Euclidean spacelike four-momenta, k'^,p'^,q^ > 0, 



they are continued to timelike four-momenta k'^ < 0, by means of an —ie prescription in the denominators of 
( 2.25[ ). Then the imaginary part of the chiral amplitude (2.7), corresponding to the cut diagram illustrated 
in Fig. [2]is given by (2.30) evaluated at s = —k'^ > 0, i.e. 



(2p -qpi+p^ p2 + q^ Pb) 



-m^po 



(2.32) 



which shows that p. = pi drops out of (2.30) for s = —k on shell, and the finite imaginary part of the 



amplitude is completely non-anomalous for timelike k'^. 

The anomaly in Re A comes about because of the cancellation of the k'^ -\- s in the denominator of the 



unsubtracted dispersion integrals (2.25) and the same factor in the spectral function sum, (2.30), resulting 
in the finite integral. 



ds p^{s]p^ 



^ -/O 



dx 



l-x 



dy 



ds 5{s — S) 



27r2 



(2.33) 



independent of p'^,q'^,m? > 0. Thus the anomalous divergence of the axial vector current is tied to a 



ultraviolet finite sum rule (2.33), for the associated spectral density Pj^{s) = pi{s). 



The finite sum rule (2.33 ) and relationship (2.30 ) between the spectral functions is critical to the infrared 
aspect of the axial anomaly. Using 2p ■ q = k'^ — p'^ — q"^ , and rearranging (2.30) we find 

r,2 



— {P2 - Pi) + — {p5 - Pi) + — Po ■ 
s s s 



(2.34) 



It is easy to see from the Feynman parameter representations (2.28) that the differences, p2 — pi and p^ — pi 
are positive for spacelike p^^cp'. Hence the function obeying the sum rule (2.33) is expressed in (2.34) as a 



sum of non-negative contributions for spacelike or null p^ and q^. If the limits p^^q^, 



m 



0~^ are taken 



(in any order), some of the spectral functions P2 — Pi, P2 — Pi and po develop logarithmic singularities, 
but each term on the right side of (2.34) multiplied by p^, g^, or approaches zero. Hence the spectral 



function p. vanishes pointwise for all s > in this combined limit. In order for this to be consistent with 



the sum rule (2.33), Pj^{s) must develop a 5{s) singularity at s = in this limit. It is straightforward either 
to calculate the function p^(s; q'^,p^) from the relations above for any p"^, q^ ^m? and verify the appearance 
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of a more and more sharply peaked spectral density in the limits, p^^q^^m? 0"'', or alternatively, to 



evaluate = pi directly from (2.28a) in this limit, where from (2.26) the function 5'(x,y;0,0) vanishes 



identically. Then by interchanging the limits and integrations over from (2.28a) and (2.31a) we obtain 

^2 



lim 



dx 







l-x 



dy 6{s) 







27r2 



6{s). 



(2.35) 



Hence the appearance of the 5{s) in this limit is explicit in this representation. This delta function 
shows that a massless pseudoscalar appears in the intermediate state of the cut triangle amplitude when 



= = Tin? = 0. 



To examine this infrared behavior in more detail, it is instructive to consider the case of 



0, 



while still retaining m as an infrared regulator. In this case from (2.26), S = m /xy^ and we easily find 



/ 2 2 




V J 



m 



po(s;0,0) 



9{s - ^-rr?) 



(2.36) 



As expected for vr? the spectral function (2.36) vanishes pointwise for all s > 0. However by making 
the change of variables s = 4m? /{I — u^) it is straightforward to verify that the integral over s of the 



function (2.36) is independent of m and given by (|2.33l). Because of \2.22\ and (|2.30l) with 2p-q = k'^, the 



spectral integral (2.25) for the full amplitude /i is 

^ r „2 
A;2 



fiik% 



e 
2t^ 



m 



ds 



k^ + s 



^(fc2;0,0;m2) 



(2.37) 



due to the sum rule (2.33). This shows that the amplitude /i develops a pole at A;^ = q when p^ = q^ 



m 



0, corresponding to the d{s), (2.35) in the imaginary part in the same limit [9j. 



When the fermion mass is non-zero the amplitude (2.37) can also be written in the form. 



2=0 



m 



(h po{s;0,0) 
4^2 s k"^ + s 



(2.38) 



which shows that the amplitude has no pole divergence as /c^ q with m? > fixed [27 . Because of 
the sum rule and relations (2.36), the residue A{k'^;0,0;m?) of the pole vanishes in this limit. This may 
be understood as a consequence of decoupling, for with no other scales remaining, the limit k'^/m? 
is equivalent to the limit m? — > cxd with k'^ fixed, in which case the entire fluctuation represented 



^ The conjectured form of the spectral function of ref. [TO], eq. (11.50), disagrees with the exact resuh, (2.36l, although the 
qualitative conclusions are unchanged. 
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by the triangle diagram should vanish on physical grounds. The spectral function representation and 
determination of the anomaly by its finite parts builds in this decoupling limit m? oo automatically. 



Conversely, if m = then the amplitude (2.37) behaves like k~ for all k , in both the infrared and 
ultraviolet, as would be expected for a conformal theory with no intrinsic mass or momentum scale. This 
decoupling behavior is also inherent in the Pauli-Villars regularization of the triangle amplitude, since the 



first anomalous first term in (2.23) is exactly the negative of the second term in the limit of infinite mass. 



For comparison we may consider /02(s), given by 
/52(s;0,0) = p5(s;0,0) 



^ 1 
27r2 s 



1 



Am? 



Am 



(2.39) 



Because of its slower fall off with s, the integral of P2{s) over s does not converge, and does not obey a 
finite sum rule. Its corresponding amplitude, 

o2 




Am^ , 
2 + + 



1 + ^ + 1 



1 + 



Am? 



1 



27r2p 



(2.40) 



does not possess a pole or even a finite limit as ^ or k"^ /m^ oo. Thus this limit, equivalent to 
fixed and A;^ — > c>o is purely ultraviolet in character, and cannot be interpreted in terms of an infrared 

massless state with a finite spectral weight. 

The full amplitude (2.13) for = = becomes 

^ = /i(fc^O,0)[rr^(p,g) +rr^(p,9)] +/2(fc^0,0)[Tr''(p,g) +rr^(p,9)] 
:/i(A:2;0,0) k^^v^^{p,q)+[h{k';0,0)-h{k^;0,0)] [rr^(p, g) + <?)] , (2.41) 



by use of the identities, (2.17). For vanishing p = q = 0, the tensors. 



ua/3 



uaf3 



-A'^"(p,g), 



(2.42a) 
(2.42b) 



have zero contraction with photon wave amplitudes obeying the transversality condition, p"^Aa{p) 



q^Afj{q) = 0. Hence the term involving the T2 + T5 in (2.41) drops out entirely in the full matrix ele- 
ment of ^5(0) between the vacuum and a physical two-photon state \p, q), giving simply 



{0\4mp,q) = iTf^''^{p,q)A^{p)Ap{q) 



p2— g2— g 



= ih{k^; 0, 0) k^^v^^^ip, q) A^{p)Ap{q) 
■ 2 

1^0 zvr^fc^ 



(2.43) 
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where the last hne fohows from (2.37) for m = 0. This exhibits the pole at k = {p + q) = 0. Thus the 
singular infrared behavior required by the anomaly, survives in the full on shell matrix element to physical 
transverse photons. The residue of the pole is determined by the anomalous divergence, 

{^\d^J^mp,q) = ik^{^\4mP:l) = v''^{p,q)A^{p)Ap{q) (2.44) 

to be non-vanishing when m = 0. 

By examining the expressions above one can see that the full amplitude exhibits propagating pole-like 
behavior for fc^ ^ p^,q^,m?, while for finite p'^,q'^,m'^, the pole appears to soften into a resonance and 
there is no singularity when k'^ < min (Ip^lJgi^l^ j^j^s-j^ Thus a strict infrared pole at A;^ = q exists in 
the triangle amplitude only for zero mass fermions, and it couples to the physical amplitude only when 



p^ = q^ = m'^ = 0, giving the full answer for the on-shell matrix element (|2.43l) of the chiral current to 



physical transverse photons only in this case [23]. However, because of the sum rule (2.33), the pseudoscalar 



state implied by the anomaly is present at any momentum or mass scale, while from (2.37) its coupling to 



(Jg) and photons simply becomes weaker for k"^ < min |^^|, w-^), and decouples entirely as fc^ ^ 
for any of \q'^\,m'^) finite. 



This appearance of a massless pseudoscalar in the chiral amplitude (2.7) in the two-fermion intermediate 
state in the limit of massless fermions is reminiscient of the Schwinger model, i.e. massless electrodynamics 
in 1 -|- 1 dimensions, where it is also related to the anomaly [6]. In each case one can use the fermion mass 
as an infrared regulator to examine the appearance of the anomaly pole in the amplitude or delta-function 
in its imaginary part as the limit tit? — > 0, for k^ < timelike. In each case when one finally arrives at 
the limit of null four-momenta, the intermediate state which gives rise to the pole is a massless electron- 
positron pair moving exactly collinearly at the speed of light [24 . Thus even the 3-1-1 dimensional case 
becomes effectively 1-1-1 dimensional in this limit, which accounts for the infrared enhancement. The only 
essential difference between d = 2 and d = 4 dimensions appears to be the necessity of going to a more 
complicated three-point amplitude in the d = 4 case to reveal the anomaly pole. The special role of the 
kinematics of the triangle diagram for this infrared enhancement in gauge theories in 3 -|- 1 dimensions has 
been emphasized previously in [r4j . 



III. THE AUXILIARY FIELD DESCRIPTION OF THE AXIAL ANOMALY AND 
ANOMALOUS CURRENT COMMUTATORS 



The appearance of a massless pseudoscalar pole in the triangle anomaly amplitude suggests that this 
can be described as the propagator of a pseudoscalar field which couples to the axial current. Indeed it 
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is not difficult to find the field description of the pole. To do so let us note first that the axial current 
expectation value ( Jg*) ^ can be obtained from an extended action principle in which we introduce an axial 
vector field, B/^i into the Dirac Lagrangian, 

i'tp-f^ (^d^i -leA^^ - mi) iip-f^ (^dfi -ieA^ - igj^B^^ ip - mi!ip (3.1) 

so that the variation of the corresponding action with respect to Bfj_ gives 

9(4) . ■ (3.2) 



Henceforth we shall set the axial vector coupling g = 1. Next let us decompose the axial vector Bf^ into its 
transverse and longitudinal parts, 

B^ = B^+ d^B (3.3) 
with d^B^ = and B a pseudoscalar. Then, by an integration by parts in the action corresponding to 



(3.1 ), we have 



W')a = -%- (3.4) 



Thus the axial anomaly (2.5) implies that there is a term in the one-loop effective action in a background 



and B^ field, linear in B of the form. 



e2 



V/ = / d^x e'^-P-F^.Fp, B , (3.5) 



or since d Bx = \3B, 

=2 



Seff = Jd'^l d^V [^^"'"F^uFp.], U-,^ [d^Bx]y , (3.6) 

where \I\~y is the Green's function for the massless scalar wave operator □ = 9^9^. Thus from (3.2), this 
non-local action gives |25] 



^"^'^^ = ^^'^□-'^"^'"^^a/^^P- , (3.7) 



which explicitly exhibts the massless scalar pole in the massless limit of (3.10), and which agrees with 
the explicit calculation of the physical (0175*1^,5) triangle amplitude to two photons (2.43) in the previous 
section for = = m? = 0. 



The non-local action (3.6) can be recast into a local form by the introduction of two pseudoscalar 
auxiliary fields t] and x satisfying the second order linear eqs. of motion, 

nr] = -d^Bx, (3.8a) 
□ X = 1^ F^.F'^-' = ^ e'^-^'-^F.^F,^ . (3.8b) 
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Then one can verify that the local quadartic action functional, 



(3.9) 



yields back the eqs. of motion (3.8) when freely varied with respect to x ^-iid t] respectively, while evaluating 



to (3.6) upon using these eqs. of motion to solve for and eliminate the auxiliary fields. In the auxiliary 



field form of the effective action, (3.9) the quantum expectation value of the chiral current (3.2) is given 
by 



4[x] 



5S, 



eff 



6B„ 



d^X 



167r2 



af3^ per 



(3.10) 



at least insofar its anomalous divergence is concerned. The effective action (3.9) reproduces the anomalous 



divergence of (J^), but not necessarily the non-anomalous parts of the tensor amplitude T^°'^{p,q). The 
r] and x fields and their propagator are a local field representation of a massless 0~ state propagating in 
the physical matrix element (2.43) with 



diagram with source FnyF^^^ on one end, and d B\ on the other, as in Fig. 3 



q = m 



0, which may be represented by the effective tree 



D 



XV 




FIG. 3: Tree Diagram of the EfFfective Action (3.9 1 showing the massless propagator D^ri representing the massless 



state in the triangle amplitude (2.431 to physical photons when the fermion mass m = 0. 



The same diagram also represents the vector current expectation value (J^)^ in the presence of a 
background axial field B\ and gauge field A^, also implied by the original triangle diagram Fig. [T] upon 
reversing the roles of the axial vertex and one of the vector vertices, i.e. 



(3.11) 



This crossing symmetry or equivalently, the fact that the non-local action (3.5) involves a mixed term 



involving both F^^F^^^ and d'^B\ is the reason why two massless pseudoscalar auxiliary fields rather 
than just one are required to describe the amplitude correctly through a local effective action. A single 
auxiliary field would necessarily produce unwanted direct F^j,uF^'^ \^~^ F^yF^^ and d'^B\\^^d^B\ terms in 
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the effective action, not present in QED. Note that because of the mixed field kinetic term in (3.9), the 
only propagator function that can appear is the mixed one, D^^ = i{Txf])- 



Several additional remarks concerning the effective action (3.9) are in order. First, the axial anomaly, 
and hence the fields r] and x ^'^^ their propagator D^rj are present for any p^.,q^-, m? , although they decouple 
from the physical amplitude T'^"^{p,q) as k"^ ^ if any one of p^,q'^,m'^ are greater than zero |10^ [26] . 
The massless rj and x fields decouple from all physical processes involving electrons in that case, and 



the amplitude has a resonant peak (or peaks) at s ~ {p ,q ,m ) as in (2.36), rather than a sharp 6{s 



behavior. Because of the sum rule (2.33) the resonance has the same total probability when integrated 



over s, but the massless propagator D^ri = iC^XV) saturates the physical on-shell amplitude (2.43), and 
may be substituted in its place only when = q'^ 



m 



0. 



Secondly, since it contains kinetic terms for the auxiliary fields r] and Xj the effective action (3.9) 



describes two massless pseudoscalar degrees of freedom. These degrees of freedom are two-particle 0~ 
correlated e^e~ states, and rj and x are pseudoscalar composite fields of bilinears of 4^ and ip. In fact, r] 
and X may be defined by their relations to the vector and axial currents and of the underlying Dirac 



theory (2.1 ) by (3.11 ) and (3.10) respectively. Hence varying rj and x aiid treating them as true degrees of 
freedom is equivalent to varying the bilinear current densities J'^ and Jg, according to (3.11) and (3.10). 



Thirdly, these pseudoscalar degrees of freedom are implied also by a canonical operator treatment. 



Taking the effective action (3.9) as defining canonical momenta conjugate to the r/ and x fields via 



5f] 

Sx 



-X: 



and imposing the equal time canonical commutation relations, 

[r/(t,f),n^(t,f)] = i53(x-f) = [x(t,x),n^(t,x')] 



we find that the currents defined by (3.10) and (3.11) satisfy the commutation relations, 

o2 



[J0(t,f),j0(t,f')] 



[JXt,f),jO(t,x')] 



e 

„2 



FO^a,[??(t,f),x(t,x)] 
,F*^9,-[r/(t,x),x(t,x)] - 



2^ 



B • V5^(x 



le 



(E X VY5^{x-^), 



(3.12a) 
(3.12b) 

(3.13) 

(3.14a) 
(3.14b) 



which are the anomalous commutation relations deducible from the covariant T* time ordering of the 
currents required by the axial anomaly [^Ql [27]. In other words, the canonical commutation relations of 



the auxiliary fields with the kinetic terms in the effective action (3.9) are required by the anomalous equal 
time commutators of the currents and J5 . This suggests that 77 and x should be treated as bona fide 
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quantum degrees of freedom in their own right. Because of the unique kinematic status of the triangle 
diagram [14' , and the non-renormahzation of the axial anomaly J28j , the structure of the effective action 



(3.9) and commutation relations (3.14) are not modified by any higher order processes. 



Fourthly, we observe that the energy corresponding to (3.9) is not positive definite. This in itself should 



not be surprising, since the action (3.9) is a finite effective action in which the formally infinite energy of 
the Dirac sea has been effectively subtracted (by the counterterms needed to impose gauge invariance, not 
encountered explicitly in our approach). Under some boundary conditions, this finite subtracted energy 
can be negative, as in the Casimir effect. The conditions under which this is true requires a careful analysis 
of the surface terms which we have neglected so far in our discussion. In fact, because both d^B^ and 



F^yF^^^ = d^{e'^'^°'^ A^Fafs), are total derivatives, the action (3.9) changes only by a surface term under 
constant shifts of either r] or Xi ^-nd there are two Noether currents, 

= ^^""""A'^F^c.p - d^x , (3.15a) 
= B^' + df'r] , (3.15b) 



with corresponding Noether charges which are conserved by the eqs. of motion (3.8). The dynamics of the 
rj and x fields are partly constrained by these conservation laws, and should be considered together with 
the dynamics of vector and axial vector sources and 



Finally, since the effective action (3.9) explicitly exhibiting these two pseudoscalar fields is nothing but 
a rewriting of the non-local form of the effective action for massless QED in the presence of an axial vector 
source, the massless degrees of freedom they represent have not been added in to the theory in an ad hoc 
manner. They are contained in QED as soon as it is extended by an arbitrary axial vector coupling as in 



(3.1), and are a necessary consequence of the axial anomaly, which in turn is required by imposition of all 
the other symmetries. 

In condensed matter physics, or electrodynamics at finite temperature or in polarizable media, where 
Lorentz invariance is broken, it is a familiar circumstance that there are low energy collective modes of 
the many-body theory, which are not part of the single particle constituent spectrum. This occurs also 
in vacuo in the two dimensional massless Schwinger model, whose anomaly and longitudinal "photon" 
can be described by the introduction of an effective scalar field composed of an e~^e~ pair [7]. In 3 -|- 1 
dimensions, relativistic kinematics and symmetries severely limit the possibilities for the appearance of 
such composite massless scalars, with the triangle anomaly the only known example [14 . The fact that 
the e~^e~ pair becomes coUinear in the massless limit shows that this effectively reduces the dimensionality 
back to 1 -|- 1. In the well studied 1-1-1 dimensional case, the commutation relations of fermion bilinear 
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currents and J5 , which create the composite e^e massless state are due to the anomaly [29j. Evidently 



from (3.13)-(3.14), a similar phenomenon occurs in the triangle amplitude in 3 + 1 dimensions. 

To conclude this section, one may ask: if there are massless pseudoscalar states in a weakly coupled 
theory like QED, which has been subjected to such exquisitely accurate tests, how could they have escaped 
detection? The answer to this is twofold. First, as we have seen these massless pseudoscalars do not couple 
to real QED with a finite electron mass, except at /c^ ^ m^, so in massive QED they have no effects on low 
energy or long range electromagnetic interactions. Second, and more importantly, they require an axial 
vector source (as well as a non-zero FF). In pure QED there is no axial vector coupling, i.e. (7 = in 



(3.1 ). Indeed it is impossible to introduce such a coupling into a U{1) gauge theory with a dynamical axial 
vector Bfj_ field without the breakdown of Ward identities necessary to the ultraviolet renormalizability of 
the theory. While it is theoretically possible to introduce a non-dynamical axial vector source, except for 
vr*^ decay where indeed the axial anomaly with a quark triangle amplitude dominates [12] , it seems to be 
difficult to realize such a source in nature, at least on macroscopic scales. In this situation the appearance 
of a massless pseudoscalar pole in the QED triangle amplitude, and its description by massless auxiliary 
fields is an interesting curiousity, illustrating the logical and kinematical possibility that anomalies may 
lead to unexpected consequences for the long distance physics in higher dimensions as well as lower ones, 
but which does not affect any predictions of QED in four dimensions with gB^ = 0. 



IV. THE (TJJ) TRIANGLE AMPLITUDE IN QED 

In this section we consider the amplitude for the trace anomaly in flat space that most closely corre- 
sponds to the triangle amplitude for the axial current anomaly reviewed in the previous section, and give a 
complete calculation of the full {T^^^ J°^jl^) amplitude for all values of the mass and the off-shell kinematic 
invariants. Although the tensor structure of this amplitude is more involved than the axial vector case, 
the kinematics is essentially the same, and the appearance of the massless pole very much analogous to 
the axial case. 

The fundamental quantity of interest is the expectation value of the energy-momentum tensor bilinear 
in the fermion fields in an external electromagnetic potential A^, 

{THa = {TZJa + (TZ). (4.1) 
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where 



(4.2a) 
(4.2b) 



are the contributions to the stress tensor of the free and interaction terms of the Dirac Lagrangian (3.1 1. 



The notations, t^^'^^ = {t^^ + t^^)/2 and = {8^i— 9p)/2, for symmetrization and ant-symmetrization 
have been used. The expectation value (T^*^)^ satisfies the partial conservation equation, 



.flU\ 



(4.3) 



upon formal use of the Dirac eq. of motion (2.2 ). Just as in the chiral case, the relation is formal because 



of the a priori ill-defined nature of the bilinear product of Dirac field operators at the same spacetime 



point in (4.2). Energy-momentum conservation in full QED (i.e. when the electromagnetic field A^^ is also 



quantized) requires adding to the fermionic T^" of (4.2) the electromagnetic Maxwell stress tensor, 

1 



Max — ^ ^ X 



HupXp p_ 



Xp 



(4.4) 



which satisfies 8uT^^j^^ = —F^^Ju- This cancels (4.3) at the operator level, so that the full stress tensor 



of QED is conserved upon using Maxwell's eqs., d^F^^ = J^. Since in our present treatment A^ is an 
arbitrary external potential, rather than a dynamical field, we consider only the fermionic parts of the 



stress tensor (4.2) whose expectation value satisfies (4.3) instead. 



At the classical level, i.e. again formally, upon use of (2.2), the trace of the fermionic stress tensor 
obeys 



™ (d) 



9pu 



T^'^^'^') = -mipi^ (classically) , 



(4.5) 



analogous to the classical relation for the axial current (2.4). From this it would appear that {T^'^)^ will 



become traceless in the massless limit m ^ 0, corresponding to the global dilation symmetry of the classical 
theory with zero mass. However, as in the case of the classical chiral symmetry, this symmetry under global 
scale transformations cannot be maintained at the quantum level, without violating the conservation law 



satisfied by a related current, in this case the partial conservation law (4.3), implied by general coordinate 



invariance. Requiring that (4.3) is preserved at the quantum level necessarily leads to a well-defined 



anomaly in the trace, [H |30l |3T], namely. 



{T> 



P' A\m=0 



24^2 f^- 



F.,.,.F^''' 



(4.6) 
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analogous to (2.5). It is the infrared consequences of this modified, anomalous trace identity and the 



appearance of massless scalar degrees of freedom for vanishing electron mass m = that we wish to study. 
Our first task is to evaluate the full amplitude at one-loop order obtained by taking two functional 



variations of the expectation value (4.1) with respect to the external potential, and then evaluating at 



vanishing external field, A = 0. In position space this is 



{ief{T^:jz)r{x)jP{y))A=o 



^''""'y^ - 6A^{x)5Ap{y) 
+6\x - z)<7"('^n-)^(z, y) + 5\y - z)/(^n'^)"(z, x) - 9^-[5\x - z) - 5\y - z)]Ii'-''{x, y) , (4.7) 



where 



™/3^ , diJ'^ix))^ 



-le 



(4.8) 



is the electromagnetic polarization tensor in zero external field. Going over to momentum space and 
factoring out the resulting factor of momentum conservation, (27r)^5^(A; — p — q), we obtain 



r^''"^(p,g)= / d^x / d^y e'f-''+''i-yT'"""^{z = 0;x,y), 



(4.9) 



which receives contributions from the two kinds of vertices derived from the two terms in (4.2), namely 



V^^ih, k2) = \ {-i^{k^ + + 7"(fci + ^2)1 - 25^"[7^(^i + ^2)a + 2m] , (4.10a) 
^/.^a = -^(7^5^" + i^g^"^) + g^^'l'' , (4.10b) 



respectively, represented in Figs. |4j 

h 





FIG. 4: The two kinds of vertices contributing to the stress tensor amplitude (4.9 1 



At the one loop level the amplitude (4.9) is represented by the diagrams in Figs. |5] together with 
those in which the photon legs are interchanged. The first of these diagrams with the vertex V^'^ gives the 
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contribution, 



V 



idH 



7 



— T 



(27r)4 [ ' ^^/+^ + m'/ + m'/-j^ + m 

2 / d^/ tr{y^'^(^+p,/-g)(-/-/5 + m)7°(-/ + m)7^(-/ + ^ + m)} 



(27r)4 [(/ + p)2 + m2] [(/ - g)2 + m2] [«2 + ^2] 

while the first of the diagrams with the vertex W^'^"' gives 



w 



idH 



tr <^ W^'"'' 



3 

— T- 



(27r)4 [ jl + rn' ^-^ + m 

2 / dH tr{VF^''°(/ + m)7^(/-^ + m)} 



(27r)^ 



[/2 + m2][(/-g)2+m2] 



where 



tr{7°(/ + m)7^(/ -|5 + m)} 
(2^ [/2+m2][(/-p)2+m2] 



(4.11) 



(4.12) 



is the Fourier transform of (4.8). The additional factor of i in the loop integration measure of (4 
(4.12) comes from the continuation to Euclidean momenta: Iq il^. 



(4.13) 
TTl) and 





FIG. 5: The two kinds of amplitudes contributing to the stress tensor amplitude (4.9 1 



As usual these loop integrals are formal and divergent, since the one-loop polarization requires regular- 
ization and renormalization, which we postpone for the moment. The second set of diagrams with the W 



interaction vertex give rise to the second set of terms in (4.7) explicitly proportional to the polarization 



(4.13 ). As we shall see, we actually require only the finite parts of the first diagram in Fig. [5| together with 
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the equivalent diagram obtained by interchanging p and q and a and [j. The contribution of the contact 
terms of the second kind of diagram in Fig. [5] with the W^^"' vertex wih be determined from the Ward 
identities. 



The fuh one-loop contribution to the amplitude (4.9) is the Bose symmetric sum 



(4.14) 



Vector current conservation (9^ J'^ = implies that both the polarization, and this amplitude should satisfy 
the Ward identities, 



p„n"^(p) = g^n"^(g) = 0; 
p«r^-°^(p,g) = g^r^'^-/3(^^^^^o_ 

The first of these relations implies that 11"^ (p) is transverse, i.e. 

n"^(p) = (p2g"/3 _ p>/3)n(p2) . 



(4.15a) 
(4.15b) 

(4.16) 



In addition, the partial conservation law (4.3) implies that the amplitude in (4.7) should also satisfy the 
Ward identity. 



^r'^-°/^(z; X, y) = n^-(y, z){g^^5t - 6^g>^^)/j6\x - z) 



+n°-(x,z)(/'^<5^ - 5P,g^^>^)^5\y - z) . 



In momentum space this becomes 

fe.r^-"/^(p,g) = (<7^V - 0'^)n^"(g) + {g'^fqu - 5Un^'"'{p) 



(4.17) 



(4.18) 



or combining with (4.16), we obtain 



q^'pyl^ - q^g^'^p^ + g^'^q'^p^ - g^^p^p • q) n{p^) 
+ fp^g"/ - p^g^'^q'^ + <?^Vg^ - 5^" A • q) n(g2) . 



(4.19) 



As already remarked, all of these relations are formal since both Ti""^'^{p,q) and n"^(p) are ill-defined a 
priori and require some procedure to extract the finite terms. 



Formally one could use (4.5) to obtain an identity for the trace on the amplitude g^^v^^'^'^^ ■ Like 



(4.15) and (4.19) this trace relation can be proven if and only if shifting the loop integration variable / 



in the integrals (4.11) and (4.12) is allowed, and terms which are odd in the loop integration variable 



are dropped. Since the loop integrals in (4.11), (4.12) and (4.13) are formally quadratically divergent 
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they are not well defined as they stand, so that formal manipulations of this kind can yield ambiguous 
or incorrect conclusions. The standard method of dealing with such ill-defined expressions is regular- 



ization. Any regularization method that preserves the Ward identities of gauge invariance (4.15) and 



coordinate covariance (4.19) may be used, such as dimensional regularization, Pauli-Villars regularization 
or the Schwinger-DeWitt heat kernel method. It is important to recognize that regularization amounts to 
supplying additional information which serves to define an ill-defined expression, by requiring that certain 
symmetries of the classical theory be strictly maintained at the quantum level. 

Here we shall follow the alternative approach, exactly parallel to the previous treatment of the axial 
triangle anomaly in Sec. 2, which does not require any explicit choice of covariant regularization scheme. 
Instead we define the Lorentz covariant tensor amplitude r^'^"^(p, g) by its finite terms, together with 



the requirement that the full amplitude satisfy the Ward identities (4.15) and (4.19). Then the joint 
requirements of: 

(i) Lorentz invariance of the vacuum, 

(ii) Bose symmetry. 



(iii) vector current conservation (4.15), 

(iv) unsubtracted dispersion relation of real and imaginary parts, and 

(v) energy-momentum tensor conservation ( |4.19 ), 

are sufficient to determine the full amplitude T^'^'^^{p,q) in terms of its explicitly finite pieces, and yield 
a well-defined finite trace anomaly. As in the axial anomaly case conisdered previously, this method of 
constructing the full r^'^"'^(p, g) may be regarded as a proof that the same finite trace anomaly must be 
obtained in any regularization scheme that respects (i)-(v) above. 

Lorentz invariance of the vacuum is assumed first by expanding the amplitude in terms of all the possible 
tensors with four indices depending on p", and the flat spacetime metric g'^'^. There are 2^ = 16 tensors 
with all four Lorentz indices {fi, v, a, 13) are carried by either p or g; 2^ x 6 = 24 tensors in which two of the 
four indices are carried by the symmetric metric tensor g^^ and the other two by either or g^; and just 
3 tensors in which the four indices are distributed over a product of two metric tensors with no factors of 



p or q. The complete set of these 43 tensor monomials is given in Table III Lorentz covariance requires 
that the amplitude r^'^°'^{p, q) must be expandable in this complete set of 43 tensors with scalar coefficient 
functions of the three invariants p^, g^, and p ■ q, or equivalently p'^,q'^, and k'^ = {p + q)'^- 



Since the amplitude (4.9) has total mass dimension 2, the scalar coefficient functions multiplying the 
tensors in our list which are homogeneous of degree 4 in p and q have mass dimension —2. These coefficients 
can be extracted in terms of loop integrals which are UV quadratically convergent and finite. Then the 
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pfJ^pt^pO^qP 


pfipi^qUqiS 


pfiqi^qOpp 


pf^qi^qaqP 






pl^p^qO^pP 


pt^qi^pO^qP 


qf^p'^qapP 


qt^pi^q^qis 






pf^q'^p'^pP 


qt^p^pO^qU 


qf^q'^p'^pP 


qt^q^p^qP 


g^>^gP^ 




qf^p'^pO^pP 






qt^qi'qO^pP 




pt^p'^gO'fi 


pPp'^gO'P- 


p0pt^gO"^ 


pO^p^gPi^ 


pf^pO^gPi^ 


pO^pPgt^^ 


pt^q^gO^P 


pPq^g^t^ 


pf>qi^gO"^ 


pO'q'^gPt^ 


pt^qO'gP'^ 


p'^q'^g^'" 


qt^p^gO^P 


qPp^g^^t^ 


qPpt^gO^^ 


qO^p'^gPt^ 


q^'p^'gP" 


q'^pfig^'" 


qt^q'^gO'P 


qPq'^gO^t^ 




qO^q'^gPt^ 


qt^qO'gP'^ 


qO^qPgt^'^ 



TABLE III: The 43 fourth rank tensor monomials into which V^^'^"!^ {p,q) can be expanded 

coefficients of the remaining tensors are determined by the Ward identities of vector current and stress- 
tensor conservation. Let us define the two-index tensors, 



w^f^ip, q) = p'^q'^g"'^ + {p ■ q)p"q'^ - ip^'p^ - p^q'^q^ , 

each of which satisfies the conditions of Bose symmetry, 

n^\p.q)=vP^{q,p), 
w'^^{p,q) = w^'^{q,p) , 

and vector current conservation. 



(4.20a) 
(4.20b) 

(4.21a) 
(4.21b) 



(4.22a) 
(4.22b) 



Pau''''{p,q)=0 = qpu''''{p,q), 
Po,w''f'{p,q)=0 = qpw''f'ip,q). 

These tensors may be obtained from the variation of local gauge invariant quantities F^yF^^ and 
{d^F\){duF''^) respectively, via 

u"^{p,q) = - 
w°'^{p, q) = 



d^x / dS (^v-^+n-y 



bA^{x)Ap{y) ' 
b''{d^F\d,F-\^)} 



bAo,{x)Ap{y) 



(4.23a) 
(4.23b) 



Making use of u°'^{p,q) and w"'^{p,q), one finds that of the 43 tensors in Table III, there are exactly 13 
linearly independent four-tensors t'^'^'^^ip, q), i = 1, . . . , 13, which satisfy 

Patr^ip, q)=0 = q^tr^P, q), i = 1, • • • , 13 . (4.24) 
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i 




1 




2 


(k^gf^" - kf'k") w°'l^{p.q) 


3 


{p^gt'" -ipt^p") u°''^{p.q) 


4 


{p^g^"'' - ipf'p") w°'>^{p.q) 


5 


(^q2gt,v _ 4^Mgi^) u°'P{p.q) 


6 


(^^2^M>^ _ 4gMq>') ty"/3(p.q) 


7 


[p ■qg'"' - 2{q^'p•' + pt^q")] u°'P{p.q) 


8 


[p ■qg'''' ~ 2(g^p'^ + p'^q^')] w'^l^{p.q) 


9 


{p ■ qp"" - p^q") [pP {qf'p" + p^^q") -p-q {gf^''p^' + /^p'')] 


10 


{p-qqP ~ q^pl^) [q" {q^'p" + pi'q'') -p-q {g^'^q'' + ff"^?")] 


11 


{p ■ qp°' - p^q"') [2 ql^q^'q'' - q^ig'^'^q'' + g'^'^q")] 


12 


[p-qq'^ ~ q^p^) [2 p°'p^'p'' - p'^{g°"'pf' + g'^^'p")] 


13 


{p^'q" + p''qf')g'^P +p- q (g"''/^ + g"^/'') - g^'^u'^'^ 
-[gP^pt^ + gPt'p'')q'^ - {g°"'q'' + g°'^'q'')pP 



TABLE IV: The 13 fourth rank tensors satisfying (4.241 



These 13 tensors are catalogued in Table IV 

This set of 13 tensors is linearly independent for generic k'^,p'^,q'^ different from zero. Five of 13 are 
Bose symmetric, namely, 



tr^P>9)=r^"(9,p), ^ = l,2,7,8,13, 
while the remaining eight tensors form four pairs related by Bose symmetry: 



(4.25) 



(4.26a) 
(4.26b) 
(4.26c) 
(4.26d) 



Expanding the amplitude (4.11) in this basis. 



13 



(4.27) 



i=l 
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Bose symmetry implies that the scalar fmictions Fi, F2, Fj, Fg, and F13 are symmetric under interchange 
of and q'^, while the remaining eight functions form four pairs related by Bose symmetry, 



Fs{k';p^q^)=F,{k';q^p'), 
F^{k^;p^q^)=F^ik^;q^p^), 
F^{e-p\q^)=Fi^{e-q\p'), 
Fn{k'-y,q')=F,2{k';q',p'), 



(4.28a) 
(4.28b) 
(4.28c) 
(4.28d) 



corresponding to (4.26). Thus there are 9 independent scalar functions in the amplitude (4.27), 5 of them 



completely symmetric, and 4 of them possessing both symmetric and anti-symmetric terms under ^ 
interchange, for 5 + (4 x 2) = 13 scalar amplitudes in all. We observe that all but ti3 contain terms which 
are homogeneous of degree four in the external momenta, whose coefficients we may constrain from the 
finite parts of the amplitude. 

We have chosen this basis so that only the first two of the thirteen tensors possess a non-zero trace. 



9,utr^P,Q)=^k'u^^ip,q), 



9,utr^{p,q)=^k'w''^{p,q), 



(4.29a) 
(4.29b) 



while the remaining eleven tensors are traceless. 



9ij.uti''°'^{p, q) = 0, i = 3, . . . , 13 . 



(4.30) 



Moreover because of eqs. (4.6) and (4.23), we have chosen the basis tj in anticipation of the result that in 



the limit of zero fermion mass, the entire trace anomaly will reside only in the first amplitude function, 
Fi{k^;p\q^). 

To proceed, we now fix as many of the 13 scalar functions Fi as possible by examining the finite terms 



in the formal expressions (4.11) and (4.12). To this end we perform the indicated Dirac algebra and 



introduce the Feynman parameterization (A2) of the product of propagator denominators. Then we make 



any necessary shifts in the loop integration variable / in (4.11) and (4.12) in order to extract only those 



terms for which the four indices (/i, u, a, (3) are carried by the external momentum vectors p and q in various 
combinations, which therefore can be removed from the loop integration. The remaining loop integration 
is then finite for these terms and can be extracted unambiguously. The details of this computation are 
given in Appendix A. 

To make this procedure completely rigorous, one can calculate first only the discontinuities of the 
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amplitude, continued to timelike four-momenta, using the Cutkovsky rule replacement, 

for the lines cut as in Fig. [2] in the chiral case. Owing to these delta functions in the cut diagram. 



the discontinuity of (4.11) is completely finite. Then for those terms (and only those terms) multiplying 
tensors of degree 4 in p and g, the real parts may be constructed from the discontinuous imaginary parts by 
unsubtracted dispersion relations, and are completely finite as well. Since there are 12 such tensors of degree 
4 in the external momenta, listed in Table [v] there are 12 finite scalar coefficient functions Cj(/c^;p^, g^) 
multiplying them which are defined in this way. The explicit form of their corresponding imaginary parts 



Pj is given by (6.2) in Sec. 6. It is not difficult to show from the linear independence of the tensors in 
Table [V| and general analyticity properties of the amplitude that the finite coefficient functions Cj of mass 
dimension —2 obtained in this way from their imaginary parts are identical to those obtained by the recipe 
of shifting the loop integration variable in the original full amplitude, in Feynman parameterized form, 
and identifying the terms multiplying each tensor listed in Table |Vj This is of course also the same result 
for the finite terms that is obtained if the loop integration were regularized in a covariant way, such as in 
the dimensional regularization or Pauli-Villars schemes, in which the shift of the loop integration variable 



is permitted. It is also noteworthy that this procedure of extracting the finite parts of (4.11) relies only 
upon the terms involving the V^"^ vertex in the triangle diagram of Fig. [5| and not the contributions of 
the W^^"^ vertex which are proportional to the polarization tensor ( 4.13[ ), and divergent. 



The 12 scalar coefficient functions listed in Table [V| are not all independent. Owing to the tensor struc- 
ture in the table imposed by Bose symmetry and vector current conservation, two pair of the coefficients 
are trivially dependent upon one other, namely 

-{p-q)Ci=q^Ci, (4.32a) 
-{P ■ q)Ci2 = P^Ce , (4.32b) 



so that only 10 independent coefficient functions can be determined from the finite parts of (4.27). 

Inspection of Table [V| shows also that the coefficients C5 and Cs are automatically Bose symmetric, while 
the remaining ten coefficients occur in five Bose conjugate pairs, viz. (Ci, C12), (C2, Cn), (C3, Cio), (C4, Ce), 
and {CjjCg), so that for example, Ci{k'^; p'^ , q^) = Ci2{k'^ q"^ , p'^) ■ Explicit formulae for all twelve finite 
coefficient functions may be given in the Feynman parameterized form, 

Cik^;p^,q^) = — C dx C dy ^^ji^^v) ^ 

' ' 47r2 _/q Jq x{1 — x) + q^ y{l — y) +2xyp ■ q + m? '' 
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j 


Cj — coefficient of 




1 




-4x2(1 -a;)(l -2a;) 


2 




— x(l — a;)(l — 4a; + Sxy) + xy 


3 


qt^q'^p^pIS 


2x{l - 2y){l - X - y + 2xy) 


4 


pi^p^pO^qP 


-2x{l - x){l - 2x){l - 2y) 


5 


{pt'q'' + qf'p'')p"qP 


x(l-a;)(l-2y)2 + y(l-j/)(l-2a;)2 


6 


qt^q^po^qp 


-22/(1 -y)(l-2x)(l- 22/) 


7 


pf^p'^q^pP 


2xy{l - 2xf 


8 


{pt'q'' + qf'p'')q'^p^ 


-2xy{l - 2x){l - 2y) 


9 


qt^qi'q^pP 


2xy{l~2yf 


10 


pl^p^qO^qP 


22/(1 - 2a;)(l - x - y + 2xy) 


11 


{pf'q'' + qt'p'')q°'q'^ 


-y(l - - 4?/ + Sxy) + xy 


12 


qt^q'^q^qP 


-4y^l - 2y){l - y) 



TABLE V: The twelve tensors with four free indices {yLvafi) on p, q which appear in the amphtude (4.271, with finite 



scalar coefficient functions Cj{k ,p ,q ) and corresponding polynomials in the Feynman parameterized form, (4.331 



where the polynomials Ci{x, y) for i = 1, . . . , 12 are listed in Table |Vj From (4.33) with the help of this 



Table it is straightforward to verify relations (4.32) and identify the Bose conjugate pairs of coefficient 



functions by interchange of x and y. These relations are verified in Appendix B. 



Identifying the coefficients of the finite amplitudes in terms of the tensors of Table IV gives 10 relations, 
which we group into the following three sets. First we have the three relations. 



Fi + 4F3 = C-j , 
Fi + 4F5 = , 



p-q {p-qY 

which multiply only the u^^ tensor. Next we have the three relations. 



C5 + C's , 



F2 + 4F4 



p2 



F2 + 4F6 = % 



F2 + 2F8 



P ■ Q 



(4.34a) 
(4.34b) 
(4.34c) 



(4.35a) 
(4.35b) 
(4.35c) 
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multiplying only the w°'^ tensor. Finally, we have the four relations, 

p-q [p-qY 

^10 — 7 ^ + 



p-q 



Fu 

Fl2 



{p ■ qy 

C3 _ C12 
Clo _ Ci 



(4.36a) 
(4.36b) 
(4.36c) 
(4.36d) 



multiplying tensors that do not appear in either u"'^ or w°'^. The first three of the relations (4.34 ) determine 
three linear combinations of the four functions Fi, F3, F5 and F-j in terms of the finite coefficient functions 



Ci, leaving only one of these four functions to be determined. Likewise the second three relations (4.351 
determine three linear combinations of the four functions, F2,F4,Fq and Fg in terms of the Ci, leaving 



only one of these four functions to be determined. Finally the last four relations (4.36) determine the four 
functions Fg, Fiq, Fu and F12 completely in terms of the Ci, and leave only Fi^ to be determined. 

The information needed to fix the remaining three functions comes from our fifth and final requirement 



on (4.27), namely the Ward identity (4.19). The contraction, k^T^^^'^^{p,q) gives six independent three- 
tensors obeying vector current conservation, and therefore six conditions on the amplitude, T^^°^l^ {p,q). 



-p^F3 + (3g2 + 4p ■ q)F^ + (2p2 + p . q)Fj - p^q^F^ - p\p^ + P ■ q)Fg + p^q^Fn = , (4.37a) 

p^Fi - (3g2 + 4p ■ q)F6 - (2/ +p.q)F8-p- qF^ + {q^ + 2p ■ q)Fn = , (4.37b) 

-p-q{p^+P- q)F9 - q^{q^ +P- q)Fii + F13 + n(/) = Q , (4.37c) 

and their Bose symmetry conjugates under interchange of p and q, 

(3p2 + 4p ■ q)Fs - q^F5 + {2q^ + p ■ q)Fy - p^q^Fg - q\q^ + p ■ g)Fio + p^q^Fu = , (4.38a) 

-(3p2 + 4p ■ q)Fi + q^FQ - {2q^ +p-q)Fs-p- qFg + {p^ + 2p ■ q)Fi2 = , (4.38b) 

-p-q{q^+p- q)Fio - p^{p^ + p ■ q)Fi2 + F13 + U{q^) = . (4.38c) 



It is evident that the symmetrized sum of (4.37a) and (4.38a) provides one new relation between ^3,^5 
and -F7, needed together with the first three relations of (4.34) to determine ^1,^3,^5 and F7 completely 
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in terms of the Cj . In this way we find 

^2 



Fi 

F5 
F7 



3 3A:2 
2C7-Cs-Cg 



12 



12A;2 



-Ci + C3 + Cs — Cg) + ( — C7 + 6*8 + Cio — C12) , 
(Ci — C3 — Cs + Cg) + Y^p- (C7 — Cs — Clo + C12) , 



-C7-C8 + 2Cb^ p 



12 

-C7 + 2C8 — Cg 



12A:2 
^2 



(Ci — C3 — Cg + Cg) + — — ^ (C7 — Cs — Clo + C 



12fc2 

,2 



12 J 



6 



+ (Ci — C3 — Cs + Cg) + ^Tp- (C7 — Cs — Cio + C12) 



+ ^C,+ 



p2C2 + g2a 



11 



{p ■ qf 



2{p-q) 



(4.39a) 
(4.39b) 
(4.39c) 

(4.39d) 



Likewise the symmetrized sum of (4.37b) and (4.38b) provides one new relation between F4,Fq and Fg, 



needed together with the second three relations of (4.35) to determine ^2,^4,^5 and Fg completely in 
terms of the C,. This gives 

C12 



F2 
F4 
Fe 
Fh 



Ci ^ C12 — Ci + 2C2 — 2C5 + 2Cii 

3(^2 3p2 3/^2 ' 

Ci 3Cio — C12 Ci — 2C2 + 2C5 — 2Cii + C12 
+ — + ^ 

Ci — 2C2 + 2C5 — 2Cii + C12 



12g2 12p2 
-Ci + 3C3 C12 



12g2 

C5 Ci 



Up' 



+ 



c 



12 



2p ■ q 6g2 6p2 



+ 



12/fc2 

Ci — 2C2 + 2C5 — 2Cii + C12 
6P 



(4.40a) 
(4.40b) 
(4.40c) 
(4.40d) 



With Fg, Fio, -Fii, -F12 given previously by (4.36) in terms of the finite coefficient functions Ci, the final 



function F13 is determined from the symmetrized sum of (4.37c) and (4.38c) to be 
n(p2)+n(g2) p2g2 ^ ^ p^a + q^Ce , 
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+ 



p- 



C5 + 



4p ■ q 



+ ^ (2C2 + C3 + Cio + 2Cii 



^ (2C2 + C4 + 2C5 + Cio) + J (C3 + 2C5 + Cg + 2Cii) 



(4.41) 



In this way all the coefficient functions Fi and hence the entire amplitude r^^"^(p, q) is determined from 
its finite parts Ci, and the one-loop polarization tensor 11, by enforcing the conservation Ward identities 
on the amplitude. In particular the trace terms involving Fi and F2 are determined unambiguously by 
this procedure, and we shall show in the next section that Fi contains the anomaly. Clearly if we had not 



enforced the conservation Ward identity relations (4.37) or (4.38), the trace is not determined, and could 
be required to satisfy the corresponding identities of conformal invariance in the massless limit removing 



any trace anomaly, at the price of violating the conservation identities (4.37) and (4.38). 



Because (4.37 ) and (4.38 ) potentially overdetermine the coefficients Fi, we note they give three additional 
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conditions on the finite coefficients Cj from tlieir anti-symmetric parts, 

/ (Ci - 2C2 + C3 - 2Cj + 2Cs) -2p-q{C7- Cg) - {2Cs - 2Cg + Cio - 2Cn + C^) = , (4.42a) 

-/ Ci + / C12 - 2p • g (/ Ci - C12) + / {-2C2 + C3 - Cio + 2Cii) = , (4.42b) 
p2(-2C2 + C4 + 2C75 + C7io) + V ■ q{-2C2 - C3 + Cw + 2Cn) + ^^(-Cg - 2C5 - + 2Cii) 



Ci + % C12 = 2 [n(p2) _ n((?2)] 



(4.42c) 



wliicli must be satisfied identically, for a consistent solution to exist. The method for verifying that these 
three conditions are indeed satisfied by the Ci is given in Appendix B. 
We note also that the difference, 

ds pu{s) 



(q'-p') 







+ s + s 



in the last of the relations (4.42) is finite, as required by the fact that all the Ci are finite. Here 

„2 



Pn{s) = -Imll 

TT 



1 + 



2m" 



4m2 „ / 9\ 
1 9{s- 4m^) , 



(4.43) 



(4.44) 



p^=-s 127r^ \ s / V s 
is the familiar spectral function of the one-loop photon polarization in QED, which tends to a constant for 



s S> Am?. Using the spectral representation of 11, the renormalization of Eqs. (4.41 ) may be accomplished 
(for non-zero m) by defining 



n^(/)=n(/)-n(o) 

so that 



ds pnjs) 
s p2 _|_ s 



2^2 



dx x{l — x) In 



1 + x{l -x) 



Pi 



i^i3 = (i^i3)H+n(o) 



with the logarithmically divergent n(0) removed by charge renormalization, 

1 1 



-R 



[i + n(o)]. 

Indeed the tensor multiplying the logarithmically divergent n(0) is 

t'ir^{p,q) = - I d^x I dS e'f"+'''-^ 



(4.45) 



(4.46) 



(4.47) 



(4.48) 



5Ao^[x)5Ap[y) ' 

which must be added to T^^'^°'^{p,q) in full QED. Thus unlike the chiral (J^J'^J^) amplitude considered 
previously, one renormalization of the amplitude (T^*^ J" J^) involving the stress tensor is necessary. How- 
ever since charge renormalization enters only through the tensor ti^ which is proportional to the stress 
tensor of the classical electromagnetic field, which is traceless, the trace of T'^'^°'^{p, q) which resides in the 
ti and t2 tensors is finite, and unaffected by renormalization. 
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THE TRACE ANOMALY AND SCALING VIOLATION 



Having determined completely the amplitude T^'^°^^{p, q) from its finite parts by the principles of Lorentz 
covariance, gauge invariance, and general coordinate invariance, we come now to the relations that would 
be expected from the classical conformal invariance of the theory in the massless limit, m ^ 0. 



The (TJJ) triangle diagram with the first vertex replaced by the naive classical trace (4.5) is 



A"^(p,g) = -m{ief 

2 r -i\ ^ ^'^^^ 
e m( — Ij 



tr 



dye 



ip-x+tq 



y{^i^r{x)j^{y)) 



(27r)^ \]t + ^ + m jl + m / — ^ + m / + ^ + m / + m ]t — + m 
This amplitude formally satisfies the conditions of vector current conservation, 



-T 



+ 



-T 



-7 



(5.1) 



]9„A"^(p,g) = = (?^A"^(p 



(5.2) 



if one is free to shift the loop momentum integration variable I'm (5.1 ). Although the integral is superficially 
linearly divergent, it is in fact at worst only logarithmically divergent because of the Dirac gamma matrix 
trace, and one factor of / is replaced by m. We may extract the factors of external momenta p and q, 
which multiply quadratically convergent integrals in d = 4 dimensions, and determine the finite parts in 



the same manner as the full amplitude. Then as before, we can determine (5.1 1 by Lorentz covariance, and 



vector current consevation, (5.2), in complete analogy to the case of the full amplitude (4.9) considered in 
the previous section. The evaluation is given in Appendix A. 

Since there are only 2 two-index tensors composed of p, q and the metric tensor, which satisfy the 



conservation conditions (5.2), namely the tensors u°'^ and w"^ defined by (4.20), we must have 

A°^(p, q) = Giu"^(p, q) + G2^x;"^(p, q) , 



(5.3) 



with Gi and G2 scalar functions of the three invariants k'^,p'^,q'^, and m^. By identifying the coefficients 
of the finite terms proportional to the four tensors homogeneous of degree two, viz. p^p^ , q'^q^, p'^q^ , and 
q"p^ , we obtain in Appendix A, 



G2 



e2 
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27r2 


Jo 






e2 




27r2 


P-q Jo 



dx 
1 







^ D 



dx 







(l-2x)(l-2^) 

ay ^ 



(5.4a) 
(5.4b) 



with D given by (2.20). Both of these vanish in the massless limit, and would be the expected values of 



the trace of the full amplitude (4.9), absent any anomalies 
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We now compare A"'^(p, q) of (5.3) and (5.4| to the exact trace of the fuh amphtude (4.27) 



9t..Tf'''''^{p, q) = Sk^Fi u^'^ip, q) + 3eF2 w'^^ip, , 



(5.5) 



computed in the previous section by requiring the Ward identities (4.37) and (4.38) of stress tensor con- 



servation. The functions Fi and i<2, given by (4.39a) and (4.40a) respectively are completely determined 
by that procedure. Let us consider the second term in (5.5) first. Since the tensor w°'^ corresponds to a 



local dimension six term in the effective action, c.f. (4.23b), we do not expect to contain any anomaly, i.e. 
we expect 3k'^F2 = G2, which vanishes in the limit m — > 0. To see this explicitly requires the following 
simple algebraic identity satisfied by the Cj coefficients, viz., 

Cs-Q-Cs-Ce + Cg + Cii = 0, (5.6) 

as is easily verified by direct substitution of Table [v| into (4.33). Subtracting twice this identity from 2>k'^F2 



given by (4.40a), and using also eqs. (4.32), we find 

1 



3k^F2 



p-q 



{p'Ci + q'C^ + 2p ■ qCs) . 



(5.7) 



In this form one may subsitute the Feynman parameterization of the coefficients (4.33), using the Table [V| 
to find 



3k'^F2 



dx dy{l-2x){l-2y)^- 

2Tr^p-qjQ Jq 

g2 ^2 .1 .l-x 



D 



2'k'^ p-q Jo jQ 



dy 



(l-2x)(l-2y) 



D 



since the first integral independent of m? in which the denominator D is cancelled. 



1 j-l-X 

dx I dy(l- 2x){l - 2y) = 
JO 



in fact vanishes identically. Comparing (5.8) with (5.4b) we verify that indeed 



2,k^Fo = Gi 



(5.9) 



(5.10) 



is non- anomalous, and vanishes in the m — > limit. 



Turning next to the Fi term in the full trace, (5.5), we need the following identity satisfied by the Cj 
coefficients, 



p2(Ci - 2C2 + C3 + 2Ci - 4C8) + (?2(2C76 - 4C8 + Cw - 2Cxx + C 
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-2p ■ q{C^ - 2C5 - C7 - + Cio) = . 



(5.11) 
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which is verfied in Appendix B. Adding this quantity to 3k'^Fi given by (4.39a) yields 



3k^Fi = p2(-2C2 + 2C3 + 2C4 + C7 - 2C8) + g2(2C76 - 2C8 + Cg + 2Cio - 2Cii) 
+2p ■ q{C3 - 2C5 + Cs + Cio) , 



(5.12) 



case, 



By substituting the Feynman parameterization integrals (4.33) and using again Table |V) we obtain in this 

{D - n?) 

(5.13) 



3A;2Fi = ^ 



1 r\—x 

dx I dy {1 — 4:xy) 











2 2 fi 



67r2 27r2 



/ dy 

JO 



(1-4x7/) 



D 



since 



1 pl—x 

dx / dy {1 — 4xy) 
JO 



1 
3 



(5.14) 



which unlike (5.9) does not vanish. 



Since the first term on the right side of (5.13) in which the factors of D cancel between numerator and 



denominator is both non-vanishing and independent of m, the trace is anomalous, and we have 

„2 



+ Gi 



(5.15) 



Hence the coefficient of the t^^'^^iip-q) tensor in the full amplitude may be written 



18^2^ 



1 -3m^ 



dx 







, (1 - Axy) \ 



(5.16) 



giving rise to the non-zero trace, 



=0 Qt,2 



(5.17) 



in the massless limit, which is exactly (4.6) in momentum space. 



It is clear that this non-vanishing trace is completely determined by the finite terms in the amplitude 
together with the imposition of the Ward identities of stress tensor conservation, for if we had not made 



use of (4.37), Fi would be still undetermined, and could be chosen to vanish in the massless limit. Of 



course, with this choice of Fi to satisfy the requirement of naive conformal invariance, the conservation 



identities of (4.37) or (4.38) would be violated, there would be an Einstein anomaly, and general coordinate 



invariance of the theory would be lost. This conflict between symmetries is quite analogous to the chiral 
case considered previously, where the naive Ward identity of C/c/i(l) invariance in the massless limit could 
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be maintained by adding an extra term to the amplitide, at the expense of violating the U{1) conservation 



identities (2.9). 



In the case of the anomalous non-zero trace of the energy-momenum tensor, it is the conformal invariance 
of the classical theory with massless fermions that cannot be maintained at the quantum level. In d = 4 
flat spacetime the conformal group is 0(4, 2) and its 15 generators consist of global dilations and the 4 
special conformal transformations together with the 10 generators of the Poincare group. The Noether 
dilation current. 



(5.18) 



is divergenceless and the corresponding charge J d^x is conserved if and only if T^^^ is traceless. A 
non-zero trace implies instead non-trivial scale dependence [U [30] . 

To see the effect of global scale transformations implied by the trace anomaly, we consider the trace 



(5.5 )a.i k = p + q = Q. Then, since 



w^"^(p, -P) = -p'u^^ip, -P) = P^p'g'"' - P>^) 



(5.19) 



and upon using (5.8) and (5.13), Eq. (5.5) for the full trace gives 

(l-4xy) + (l-2x)(l-2y) 



+ 



2 2 !•! 



'67r2 + 27r2 



2 2 !•! 



du 



dx / dy 
Jo 

u{l — u) 



p^[x + y){l — X — y) + m? 

(pV^-pV), 



(5.20) 



67r2 ' vr^ _/q p^u{1 — u) + m? 
after changing variables io u = x + y and v = {x — y) /2, and integrating over v. At the same time we note 



from (4.45) that 



2p 



dp^ 



67r2 



+ 



2 2 fl 



du 



x{\ — x) 



lo p'^x{\ — x) + m? 

is related to the /3 function of the electromagnetic coupling e^. Comparing with (5.20), we secure [8] 



(5.21) 



5..r'^^"^(p, -P) = - py') = -^^(/^"^ - pV) 



(5.22) 



a result that remains valid to all orders in perturbation theory [30]. Eq. (5.22) may also be derived by 



differentiating the Ward identity (4.19) with respect to (or p^), and then setting q = —p. Hence the 



breaking of scale invariance by the trace of the energy- momentum tensor, together with its conservation, 
may be regarded as responsible for the f3 function running of the coupling, without any direct reference to 
ultraviolet renormalization. 
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We note also that by combining the two terms in (5.21), the one-loop (3 function here, 

-1 



^ 2 



dx 



X (1 - xY 



(5.23) 







p^x{\ — x) + m? 307r^ m? ' 

vanishes as — > 0. This is the correct behavior for the infrared running with of the physical renormalized 
coupling at momenta small compared to the electron mass, where vacuum polarization is negligible and 
decoupling of the electron loop must hold [32]. For ^ m^, (3{e'^;p'^) — > e^/Gvr^, which is then identical 
to the usual ultraviolet (3 function, calculated e.g. in dimensional regularization where the infrared mass 
plays no role. We see then the physical necessity of the trace anomaly in a different way, for if the constant 



first term in (5.21) determined by the trace anomaly at m = were not present, there would be nothing 
to cancel the second integral as p"^ — > 0, and decoupling of heavy degrees of freedom at large distances 



(p <C m ) in (5.23) would not occur. 



VI. SPECTRAL REPRESENTATION, SUM RULE AND THE MASSLESS SCALAR POLE 



The physical meaning of the anomaly is further exposed by considering the spectral representation of 
the amplitude, cut across two of its legs as in Fig. [2] Following the pattern of the chiral case considered 



previously in (2.28), the spectral representations for the amplitudes Cj may be introduced by using the 

Pj{s;p'^,q'^) 



definition (2.26) and the identity (2.27), to obtain 



Cj{k';p',q') 



ds 



k^ + s 



(6.1) 



with 



Pj{s;p'^,q'^) 



47r2 



\x r ^ dy 6 {s - S{x, y; p\ q^)) 



(6.2) 



/o JO 

The Pi defined in this way are not necessarily positive, nor are they independent, owing to the relations 



(4.42 ), (5.6 ), and (5.11 ). Indeed the pj satisfy exactly the same identities as the corresponding Cj, of which 
they are just the discontinuity or imaginary part as k'^ is analytically continued to —s with s > 0. Of 



interest to us however is only the linear combination which appears in the trace (5.12). By using (5.11), 
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and repeating the steps that led from (4.39a) to (5.13), we obtain with the help of (2.26) and table [V] 
{P7 + P& + pg) + / {-Pi + Pa + /08 - pg) + {-pi + Ps + Pio - Pu) 



: p^(-2p2 + 2p3 + 2/)4 + P7- 2p8) + q^{2pQ - 2p8 + pg + 2pio - 2pu) + 2p ■ q{p3 - 2p5 + ps + Pio) 



dx 
^0 



xy 



{p\-2c2 + 2C3 + 2C4 + C7 - 2C8) 



+ g2(2c6 - 2C8 + eg + 2cio - 2cii) + 2p • g(c3 - 2C5 + C8 + Clo)} 

(^(s - S) 



1-x 



2/ '^^ 

27r^ Jo JO 
2n^ ■ 







xy 

dy 5{s - S){1 - Axy) 



{l-Axy) (D-m^) 







xy 



{k + s)p^ - m pr, 



(6.3) 



where 



/)j,(s;p ,q ) = Pa - 2p5 + P8 + P10 

,2 /•! rl-x 

dx I dy {1 - 4xy) 6 {s - S{x,y;p'^,q'^)) 



2^ 



and 



/ 2 2n 

Pm(s;p ,q 



JO 
„2 /■! 



27r^ ./o JO 



1 (1 - 4X7/) . / 2 2nN 

dy d[s- S{x,y;p ,q )) , 



(6.4a) 
(6.4b) 



since C3 — 2c5 + C8 + cio = 2xy(l — 4xy). 



These relations may be compared to their somewhat simpler analogs, (2.30) and (2.31), in the chiral 



case. Since 1 — 4xy > over the indicated range of x, y, both p^. and pm are non-negative functions of s for 



spacelike p'^ and q'^. Notice that at k'^ = —s the quantity drops out of (6.3), so that the discontinuity 



or imaginary part of (5.13) vanishes in the conformal limit m — > 0, and is non- anomalous. 



As in the chiral case (2.33), we find that spectral function which determines the anomaly satisfies a sum 



rule 



ds Pj,{s;p'^,q'^) 



-1 j-l—x POO g2 

, dx / dy(l-4xy) / ds 5{s - S) = 
2vr^ ./o 7o Jo 67r^ 



(6.5) 



by (5.14), which is independent of > 0,g^ > and m? > 0, since then S{x,y;p'^,q'^) > and the 5 
function can be satisfied over the range of s > 0. On the other hand, using 2p ■ q = k'^ — p'^ — g^, and 



rearranging the second and last lines of (6.3) gives 

„2^ _ 



Pri^'^P >^ ) = — (-2p2+P3+2p4 + 2p5+P7-3p8-Pl0) 

s 



q^ TTi^ 

H (-P3 + 2p5 + 2pQ - 3p8 + P9 + Pio - 2pii) H p„ 

s s 



(6.6) 



Since the pi develop at worst logarithmic singularities in the combined limit p'^,q^, vn? — > 0+ (taken m any 



order), (6.6) shows that p^ vanishes pointwise for all s > in this limit. The only way that this can be 
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consistent with the sum rule (6.5) is if Pj, develops a 5 function singularity at s = in this limit. Indeed 



since the function S {x , y; p"^ , q'^) defined by (2.26) vanishes identically in this limit, we see directly from 
( 6.4a[ ) that 

g2 H rl-x g2 

, lim ^PT(■s;P^9^) = / dx dy {1 - 4xy) 6{s) = S{s) , (6.7) 

by taking the limits inside the integral. Thus p^, exhibits a massless scalar intermediate state in the 
two-particle cut amplitude. ^ 

It is instructive to retain the non-zero fermion mass m > as an infrared regulator, in order to study 



this intermediate state in more detail. Comparing with (2.36) from the axial anomaly, we find when 
= q^ = that 

TOp' 1 / 4777.'^ \ 



27r2 s2 

This function is plotted in Fig. |6j 



4m^ 
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4m2 
s 


1-1 


A- 


Am? 
s 



is - ^rr?) 




(6.8) 



FIG. 6: The spectral function of (6.8 1 in units of m ^ as a function of 



4m2 ■ 



The corresponding amplitude is 

i^i(A:';0,0) 

_ 1 

~ 3P 



ds 



3fc2 Jq k^+s 



[(/c^ + s)p^ - rr?pm] 



^_^2 r Pm{s;0,0) 
67r2 Jq k'^ + s 



(6.9) 



^ This observation was made in ref. [34l in the context of photon pair creation by a cosmological gravitational field. 
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which exhibits a pole at A:^ = when the fermion mass m = 0, i.e. 



hm 



1 



187r2 ■ 



When the fermion mass is non-zero the ampUtude (6.9) can also be written in the form, 

,2 



m 



ds pm{s;0,0) 

Am? "5 



fc2 + 



(6.10) 



(6.11) 



which shows that there is no pole divergence as A;^ — > with > fixed. This is again the phenomenon 
of decoupling, as in the cancellation noted previously in (5.21) and (5.23), with the place oip^ <^ m? being 



taken here by k ^ m . Conversely, if m = the amplitude (6.9) behaves like k~ for all k , in both the 
infrared and ultraviolet, as expected in the classical conformal limit of a theory with no intrinsic mass or 
momentum scale. 

It is also instructive to carry out the same steps for the imaginary part of the trace which is non- 
anomalous, i.e. for the spectral function corresponding to the non-anomalous amplitude F2. Repeating 



the steps which led from (4.40a) to (5.8), using the identities corresponding to (4.32) and (5.6) for their 



imaginary parts, we obtain 



^2 j^2 

-2 Pi + ^ P12 - Pi + 2/02 - 2p5 + 2pu - P12 



e2 1 



p-q 



2TT^p-q Jo 



e2 1 



27r2 



P-q Jo 



dx 



dx 



^ dy ^^^Z7^ (1 - 2a;)(l - 2y) {D - m^) 







xy 



l-x 



dy 5{s - S){1 - 2x){l - 2y) 



k^ + s 



m 
xy 
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In this case the spectral function corresponding to Pj, in (6.4a) is 



p ■ q 



1 



27r2 p ■ q 



dx 











dy (1 - 2x)(l - 2y) 6{s - S{x, y;p^, q^)) 



but unlike (6.5) p% obeys the vanishing sum rule 

o2 



ds psis; p'^ , q'^^ 







e 



dx 



dyil-2x){l-2y) = 0. 



(6.12) 



(6.13) 



(6.14) 



Hence, although pg can be expressed in the form 
a2 ^2 



P8{s;p- 



2 „2\ 



p 



(P4-P8) + — (Pe-Ps) 
s s 



m? 
~2^ 



dx 



(0.15) 







xy 



analogous to (6.6), which vanishes pointwise in the combined limit, p ,q ,m 0~^, it has no positivity 
property, and no reason to develop a 6 function singularity at s = in that limit. Indeed it is not difficult 
to see that ps and indeed the corresponding full amplitude Cg vanishes identically in this limit, c.f. eq. 



(6.19b) below, consistent with the vanishing of F2, which unlike -Fi has no pole as p ,q ,m — > O"*" |35] . 
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We may use the general amplitude rf^'^°'^{p, q) to evaluate the matrix element of T^*^ to physical photons 
on shell, = = which are also transverse. In this case all terms with and vanish when contracted 
with the transverse photon polarization states, and the matrix element simplifies considerably. The tensor 
w°^^{p, q) and hence the tensors t2,t4, Iq and vanish upon contraction with transverse photons, as do the 
tensors tg, tio, in and ti2- The remaining relevant form factors also simplify considerably when = (f' = 0, 
becoming 

2Cj + Cg 



i^i(r;0,o) 



F3(/c2;0,0) = F5(A;2;0,0) 



p2—q2—Q 

Cj — Cs 



12 



i^i3,ij(A:';0,0) 



4F3(A:^;0,0) = 



-C7 + Cs 



p2— q2— Q 



2=0 ■ 



for any m. The scalar coefficients Cj here are evaluated on the photon mass shell, given by 



e 

4^ 



dx 



(6.16a) 
(6.16b) 
(6.16c) 
(6.16d) 

(6.17) 



Q xy + m? ' 

in which case C2 = Cn, C3 = Cio, and C7 = Cg. Thus, for on-shell photons there are only three 
independent form factors, and we can write the matrix element to physical photons with transverse field 
amplitudes ^a(p), Ap[q) in the form, 

(0|r^''(0)|p,g) = F, {k'g'^'' - k^^k^) u^" {p,q)A^{p)Ap{q) 

-2F3 [k" - A^p^'q^ + q'^p^) + 2(pV + q^ql] «"^(p, q)A^ip)Ap{q) 
+ F^^Rt>tr^{p,q)A^{p)Ap{q), 



(6.18) 



with Fi, F3 and Fi^r evaluated at = = given by (6.16), (6.17) and Table [v} The survival of only 3 
independent tensors when both photons are on their mass shell and have physical transverse polarizations 



agrees with the literature |3H l35] . Each of the three tensors remaining in (6.18) is conserved and their 
contractions with k^ vanish for photons on shell. Only the first has non-zero trace. 
Taking the m = limit gives the further simplification that 



C7(fc^)|p2=g2=„2=o 

C8(A;^)|p2=g2=^2=o 



e 

2^r2p 
0, 



dx 



l-x 



dy (1 - 2x) 



e' 1 
127r2 P 



(6.19a) 
(6.19b) 



while Fi?jR contains a \'a.{k'^ / rn?) behavior in this limit, but no pole (reflecting the need to renormalize the 
charge at a mass scale /i^ > different from m? in the massless limit). Thus from (6.16) both the -Fi and 
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form factors of the scattering amplitude to physical on shell photons exhibits a massless scalar pole in 



the limit of vanishing electron mass, with Fi given by (6.10) and 

e2 1 



lim 2F3(r;0,0) 



m— >0 



72^2 1^2 ■ 



(6.20) 



The leading order behavior as k'^ 
photons is 



of the sum of terms in the amplitude (6.18) to physical on shell 



lim lim {0\T'"'(0)\p,q) ^ (pf'p'' + qf'q'') u°'i^(p,q)AJp)A0(q) + log(k^) & finite terms, (6.21) 

when the limit of vanishing electron mass is taken first. This shows that the singular massless pole behavior 
survives in the matrix element of the stress tensor to physical transverse photons (in its tracefree terms), 
while the trace remains finite in the conformal limit of vanishing electron mass and all 3 four-momenta 
{k,p,q) becoming lightlike. 



The kinematics of the state appearing in the imaginary part and spectral function (6.7) in this limit is 



essentially 1 + 1 dimensional, and can be represented as the two-particle collinear e~^e pair in Fig. [7| This 
is the only configuration possible for one particle with four-momentum k^ converting to two particles of 
zero mass, p^ = q^ = as k^ ^ as well. A detailed examination of the imaginary part of the amplitude, 
illustrated by the analog of Fig. |2] shows that there is a cancellation between the numerator and Feynman 
propagator in the denominator of the amplitude from the uncut fermion line in the triangle. Thus all 
particles in the real propagating intermediate state depicted in Fig. [7] are massless, on shell, and collinear. 
Although this special collinear kinematics is a set of vanishing measure in the two particle phase space, the 



6{s) in the spectral function (6.7) and finiteness of the anomaly itself shows that this pair state couples to 



on shell photons on the one hand, and gravitational metric perturbations on the other hand, with finite 
amplitude. When gravitational scattering is considered in Sec. 8 the four-momentum transfer k^^ may be 



timelike or spacelike, the pole terms (6.9) and (6.20) in the real part of the amplitude become relevant. 



and neither fermion pair nor final state photons are collinear. 



J 



mm 



FIG. 7: The two particle intermediate state of a collinear e^e pair responsible for the S-fn. in (6.7 1 
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VII. ANOMALY EFFECTIVE ACTION AND MASSLESS SCALAR FIELDS 



Having demonstrated the existence of a real massless spin-0 intermediate state in the imaginary part 
of the triangle amplitude, and a corresponding massless pole in the full amplitude, we turn in this section 
to the effective action and scalar fields which describe these massless scalar degrees of freedom. In fact, a 
covariant action for the trace anomaly in a general curved space has been given already in several earlier 
works jTSl [36| 1371 138] . This general effective action can be presented in the non-local form, 

Sano'm[g,A] = (7.1) 



^ j d^x^g j d^x'^'{E-'^-UR^ ^l\x,x') 



2hC'^ + b' { E - - UR] + 2cF^^F^"' 



where the b and b' parameters are the coefficients of the Weyl tensor squared, 
RxuupR^^"" - 2R.^R^'' + ^ and the Euler density E =*Rxavp*R^^''P = RxuupR^^"^ - ^RnuR^" + R^ 



respectively of the trace anomaly in a general background curved spacetime, and the last term in (7.1| 



takes into account the anomaly in a background gauge field with coefficient c. For the present case of Dirac 
fermions in a classical gravitational {g^u) and classical electromagnetic (Aa) background, b = l/3207r^, and 
b' = — ll/57607r^, and c = — e^/247r^. The notation A^^{x, x') denotes the Green's function inverse of the 
conformally covariant fourth order differential operator defined by 

A4 = Vf, (v'^V' + 2R'"' - '^Rg^'^^ = + 2R^"'Vf,V^ + ^(V^i?)V^ - '^RO. (7.2) 



By varying (7.1) multiple times with respect to the background metric g^j^^ and/or the background gauge 
fields A a one can derive formulae for the trace anomaly related parts of amplitudes involving multiple 
insertions of the stress tensor {TTT...J J) and {TTT...) in curved or flat space. We emphasize that the 



effective action (7.1) was obtained by integrating the anomaly, and is determined up to terms which are 



conformally invariant. Therefore one can expect it to yield correct results for the trace related parts of 



amplitudes such as (4.11), while the tracefree parts are not given uniquely by (7.1) 



As detailed in ref. [T^ we may render the non-local anomaly action (7.1) into a local form, by the 



introduction of two scalar auxiliary fields ip and -0 which satisfy fourth order differential eqs.. 



■OR 



(7.3a) 
(7.3b) 



where we have added the last term in (7.3b) to take account of the background gauge field. This local 



effective action corresponding to (17. 1|) in a general curved space is given by 



Sa 



" '-'anom ^ ^'-'anora 2 / 



(7.4) 
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where 



nrt 



(fx- 



R 



3" ; v-Mr.v.r. ^ 3. 



+ 2 CA^.pC^^^"^ +Ue-'^-ur]^ 



(7.5) 



The free variation of the local action (7.4)-(7.5) with respect to tj^ and ip yields the eqs. of motion (7.3). 



Each of these terms when varied with respect to the background metric gives a stress-energy tensor in 



terms of the auxiliary fields satisfying eqs. (7.3). Since we are interested here in only the first variation 



of the action with respect to g^u, we may drop all terms in (7.4) which are second order or higher in the 
metric deviations from flat space. Also, if we solve (7.3b) formally for ijj in flat space, we find a D"^ k~'^ 



pole in this stress tensor. The simplest way to eliminate this higher order pole is to assume that <p is 



also first order in metric deviations from flat space, so that the entire b'S^^^ contribution to (7.4) can be 



neglected as well. These reductions are equivalent to replacing the general non-local effective action of the 



anomaly (7.1) by the much simpler form 



Sa 



I <fx^g j d^x'^'R^ , 



valid to first order in metric variations around flat space, or its local equivalent ( 7.4 ) by 

R 



d^x 



-9 



(7.6) 



(7.7) 



where 



c 
2 



(7.8a) 
(7.8b) 

□ = (7.8c) 
Then after variation we may set 93 = in flat space, and the only terms which remain in the stress tensor 



R 
"3 



derived from (7.4) are those linear in ip' , viz. 



-9^9tiu{z) 



(7.9) 



flat,ip=0 

which is independent of b and b' , and contain only second order differential operators, after the definition 



(7.8a). Solving (7.8b) formally for -0' and substituting in (|7.9|), we find 

c 



anom \^ I 



(7.10) 
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a result that may be derived directly from (7.6) as well 



By varying (7.10) again with respect to the background gauge potentials, making use of (4.23a) and 



Fourier transforming, we obtain 

^Z:^{P.a) = I d^x I d^ye^P-^^^-y = ^l^p {a^^k^ - k^k^) n-^P^l) , (7.11) 

which coincides with the first term of (4.11), with ( 6.10| ), and gives the full trace for massless fermions, 

9,uT^nom = cF^pF'^^ = -^F^pF^P , (7.12) 



in agreement with (4.6). We observe that as in the chiral case, the strict 1/k pole in the anomalous 



amplitude Fi obtained from the D^i^ = i{T%l:'ip) propagator auxiliary field applies only in the limit of 
(6.10), or equivalently for k"^ S> \q^\.,m?). The spectral representations and sum rule of the previous 



section show that when this condition is not satisfied, the two-particle intermediate state in the anomalous 
amplitude becomes a broad resonance instead of an isolated pole, as in Fig. [6j and the residue of the pole 



at A;^ = vanishes when any of p^, or are non-zero. The tree amplitude of the effective action (7.7) 
which reproduces the pole in the trace part of the {TJJ) triangle amplitude is illustrated in Fig. Isl 




FIG. 8: Tree Diagram of the Effective Action (7.7 1, which reproduces the trace of the triangle anomaly. The dashed 



line denotes the propagator D^iy, = □ ^ of the scalar intermediate state, while as in Fig. [Tjthe jagged line denotes 
the gravitational metric field variation h^^ = Sg^^. 



Most of the remarks about the auxiliary field description of the axial anomaly at the end of Sec. 3 



apply also to the trace anomaly case. As in the case of the axial anomaly, the effective action (7.5 ) or (7.7) 



explicitly exhibiting these two scalar fields is a rewriting of a part of the non-local form of the effective 



action (7.1) or (7.6) for massless QED in curved spacetime, with the reduction to (7.6) correct to leading 
order in the metric deviation from flat space, h^i, = Sg^^. The massless degrees of freedom ip and ip' are 
a necessary consequence of the trace anomaly, required by imposition of all the other symmetries. In this 
case these are scalar rather than pseudoscalar degrees of freedom. As in the chiral case, two independent 
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fields are required, and tlie propagator appearing in the intermediate state of the triangle amplitude is a 



certain off-diagonal D^i^ = i{Tip'ip) term. Unlike the chiral case the general effective action (7.5) or (7.7 1 



requires the fourth order differential operator A4 of (7.2), implying that higher order amplitudes such as 
(TTJJ) should have double poles. 

An important physical difference with the axial case is that the introduction of a chiral current J5 and 
axial vector source corresponding to it appear rather artificial, and difficult to realize in nature, whereas 
the trace of the stress tensor obtained by a conformal variation of the effective action is simply a particular 
metric variation already present in the QED Lagrangian in curved space, required by general coordinate 
invariance and the Equivalence Principle, without any additional couplings or extraneous fields. Since the 
stress-energy tensor couples to the universal force of gravity, we should expect that physical processes can 
excite the scalar 93 and ip' scalar degrees of freedom required by the trace anomaly with a gravitational 
coupling strength. If m = these produce effects of arbitrarily long range. An example of this coupling to 
a gravitational scattering amplitude is given in the next section. 



Finally we remark that strictly speaking, the anomaly action (7.1) or (7.7) and stress tensor derived 
from it contain no information about the non-anomalous or tracefree amplitudes Fi,i = 3, ... ,13, although 



in certain circumstances the addition of homogeneous solutions to the wave eqs. (7.3) can give tracefree 
terms in the stress tensor which have physical consequences |18] . Our detailed computation of the full 



amplitude (4.11) shows that there is also a massless pole appearing in the traceless part of the physical 



amplitude to two photons, (6.18) and (6.20). This traceless pole term corresponds to a term in the effective 
action of the form, 



6 



3 {df'd''F^(i)F'''^ + ] (5^^ 0-4 a'^a^) F^pF 



(7.13) 



The tensor structure of this term precludes writing it as a scalar particle exchange. The pole in this ampli- 
tude with non-trivial tensor structure is clearly connected with the possible non-zero values of d^d^F^p in 
the background electromagnetic field, which breaks Lorentz invariance. Thus it appears that in this case of 
a non- vanishing background field which is non-gravitational in origin, the Ward identities obeyed by the full 



amplitude (4.11) implies the existence of additional massless intermediate states which do not transform 
as spacetime scalars, and therefore cannot in general be described by the anomaly induced effective action 



(7.5) or (7.4). Instead these massless modes are associated with longitudinal components of the metric 
perturbation in a Lorentz non-invariant background, analogous to longitudinal plasmon excitations in a 
finite temperature electromagnetic plasma. 
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VIII. SCALAR ANOMALY POLE CONTRIBUTION TO GRAVITATIONAL SCATTERING 

In order to verify the existence of the massless scalar pole in a physical process, we consider the 
simple tree diagram of gravitational exchange between an arbitrary conserved stress-energy source T''^'^ 
and photons illustrated in Fig. l9j 




FIG. 9: Tree Level Gravitational Scattering Amplitude 



This process is described by the scattering amplitude |39] 



M = 8irG / d'^x' / d'^x 



rf'^x') 



U L, 



r^(x') 



nix) 



(8.1) 



The relative factor of — ^ between the two terms is dictated by the requirement that there be no scalar 
or ghost state exchanged between the two sources, and is exactly the prediction of General Relativity, 
linearized about flat space. That only a spin-2 propagating degree of freedom is exchanged between the 
two sources in Fig. [9] can be verified by introducing the following 3 + 1 decomposition for each of the 
conserved stress tensors. 



nOO 



Too , 



1 

V2 



<9* ^ Too : 



rj^ij ^ rj.±ij ^ Qi^ y±j ^ gj ^ y±i 



1 



+ 



B.2a) 
^.2b) 

3.2c) 



where diV-^^ = 0, diT-^^^ = T^ = 0, and V"^ denotes the static Gr een's function of the Laplacian 
operator, = d^di in flat space. This parameterization assumes only the conservation of the stress-tensor 
source(s), i.e. d^T^^'^ = 0, so that there remain six independent components of T^'^ which must be specified. 
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and we have chosen these six to be Too , T^"*" * , T-*- ^-^ and the total trace T'^, which is a spacetime scalar. 



Substituting the decomposition (8.2) into (8.1) gives 



M = SttG I (fx' I (fx 
1 



T' 
-in 



00 



(V 



2^2 



□ ^00+2 ^00 



2V: 
1 



1 

/x' 



X' ,x 

1 



I T'^ -I- - T'^ 



1 

V2 



00 



which becomes 



1 



2K- 



1 



1 



•4) 



in momentum space. These expressions show that only the spatially transverse and tracefree components 
of the stress tensor, exchange a physical propagating helicity ±2 graviton in the intermediate state, 
characterized by a Feynman (or for classical interactions, a retarded) massless propagator — D"^ — > k~'^ pole 



in the first term of (8.3) or (8.4). All the other terms in either expression contain only an instantaneous 
Coulomb-like interaction — V"^ — > k~'^ or — > k~'^ between the sources, in which no propagating 
physical particle appears in the intermediate state of the cut diagram. This is the gravitational analog of 
the decomposition, 

J' = P, 

r = j^' 



1 



of the conserved electromagnetic current and corresponding tree level scattering amplitude. 



B.5a) 
^.5b) 



d^x' I d'^xJ'^ix') { ^ ) J/,(x) 



1 



A;2 



k^ 



i.6) 



which shows that only a helicity ±1 photon is exchanged between the transverse components of the current. 



the last term in (8.6) being the instantaneous Coulomb interaction between the charge densities. 



We now replace one of the stress tensor sources by the matrix element (6.18) of the one-loop anomalous 
amplitude, considering first the trace term with the anomaly pole in Fi. This corresponds to the diagram 
in Fig. [TOj We find for this term, 

{0\Too\p,q), = -J?F^{k')u''^{p,q)A^{p)Ap{q) 
{^\K\p,q), = 3k^F^ik^)u''^{p,q)A^{p)Ap{q) 
{0\Vr^\p,q), = {0\T,j\p,q),=0. 



B.7a) 
^.7b) 
3.7c) 



Hence the scattering amplitude (8.4) becomes simply. 



Ml = ^^GT\Fi{k')u''P{p,q)A^{p)Ap{q) 



AttG 



3 "k^ 



.8) 
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FIG. 10: Gravitational Scattering of Photons from the source T'^"^ via the triangle amplitude 



where (6.9) has been used. Thus for massless fermions the pole in the anomaly amplitude becomes a scalar 



pole in the gravitational scattering amplitude, appearing in the intermediate state as a massless scalar 
exchange between the traces of the energy-momentum tensors on each side. The standard gravitational 
interaction with the source has produced an effective interaction between the scalar auxiliary field ij:' and 
the trace T'^n with a well defined gravitational coupling. Thus we may equally well represent the scattering 



as Fig. 10 involving the fermion triangle, or as the tree level diagram Fig. 11 of the effective theory, with 
a scalar particle exchange. 





FIG. 11: Gravitational scattering of photons from the trace of a source T'^ via massless scalar exchange in the 



effective theory of (8.9) 



This diagram is generated by the effective action in flat space modified from (7.7) to 

d'^x \f--g 



.9) 
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to include the coupling to the trace of the energy-momentum tensor of any matter T'^^ source. Corre- 



spondingly the eq. (7.8c) for 99 becomes 



8vrG 



.10) 



instead (7.8c). The eq. of motion for remains (7.8b) and is unaffected. We note that if the source T'^'^ 



generates the curvature R by Einstein's eqs., then R = —SttGT''^^, so that (8.9) and (8.10) are equivalent 



to (7.7) and (7.8c) at leading order in G. 



We conclude that in the conformal limit of massless electrons, the pole in the trace sector of the anomaly 
amplitude contributes to gravitational scattering amplitudes as would a scalar field coupled to the trace 
of the energy-momentum tensor of classical sources. The gravitationally coupled intermediate scalar can 
be understood as arising from collinear e~^e~ correlated pairs in a total spin 0^ state. Although the result 
appears similar in some respects to a Jordan-Brans-Dicke scalar [30] , the coupling induced by the anomaly 
involves two scalar fields each coupling to a different source, with an off-diagonal propagator, D^'^. Hence 
the phenomenology of this scalar coupling will be quite different, and the observational limits on a Jordan- 
Brans-Dicke scalar do not apply jJT] . In particular there is no direct coupling of classical energy-momentum 
sources T'|^ to T"j^ via scalar exchange as there would be in a classical scalar-tensor theory. 

Another important difference is that as we have seen, the anomaly pole is a necessary consequence of 
quantum fluctuations and low energy symmetries, whereas in classical scalar-tensor theories a postulated 
scalar field is simply added to Einstein's General Relativity. As a consequence there are one or more 



free parameters introduced in such an approach, whereas the effective action (8.9) is completely specified 



without any arbitrariness or free parameters, once the underlying quantum theory's matter content and 
couplings are given. It will be interesting to study the consequences for astrophysics and cosmology of this 
effective action derived from quantum first principles and fundamental low energy symmetries. 

IX. SUMMARY 



We have presented a complete calculation of the {TJJ) triangle amplitude in QED, for all values of the 
kinematic invariants and electron mass. As a consequence of the trace anomaly, this amplitude exhibits 
a massless pole in the conformal limit, which contributes to long range gravitational interactions, and is 
associated with the exchange of a massless 0"*" degree of freedom. This scalar exchange is described by a 



low energy local effective action (8.9) with two massless dynamical scalar fields (p and -0'. 



For the benefit of the reader we provide here a summary of the main results to be found in each section 
of the paper. 
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We reviewed in Sec. 2 the derivation of the axial anomaly in QED, showing how the finite parts of 
the triangle amplitude, together with the symmetry principles of Lorentz invariance, gauge invariance, and 



Bose exchange symmetry are sufficient to yield the complete amplitude, (2.7) for any mass and any value 
of the kinematic invariants, and determine the axial anomaly, without any explicit need of regularization 
of ultraviolet divergent integrals. We showed that the anomaly is closely connected to a finite sum rule 



of the spectral density (2.33) obtained by cutting the amplitude as in Fig. |2j For physical, transverse 
photons on shell, this spectral density vanishes pointwise for all s > 0, becoming proportional to 5{s) in 



the conformal limit of massless fermions, (2.35). Corresponding to the 6{s) in the spectral weight is a 



massless pseudoscalar pole singularity in the full amplitude and matrix element to physical photons (2.43). 
This illustrates the infrared aspect of the anomaly, and the appearance of massless states in a theory with 
anomalies in 3 + 1 dimensions. 

We showed next in Sec. 3 that the anomaly pole in the chiral case implies the existence of a non-local 



effective action, (3.6) which can be brought into a local form by the introduction of two pseudoscalar 



auxiliary fields (3.8)-(3.9). These fields and the anomaly pole can be understood as arising from a certain 
correlated two-particle collinear e^e~ state in the massless limit. The local effective action of the auxiliary 
fields has kinetic terms, and their canonical commutation relations reproduces the anomalous current 



commutation relations of the underlying fermionic theory (3.14). Thus the auxiliary fields appear to be 
bona fide massless pseudoscalar degrees of freedom, required by the chiral anomaly. 

In Sec. 4 we presented a full computation of the {TJJ) triangle amplitude in QED, where the chiral 
current is replaced by the fermionic energy-momentum-stress tensor T^'^. Following the same method as in 
the chiral case, we showed how the finite parts of the triangle amplitude, together with the same symmetry 
principles of Lorentz invariance, gauge invariance, and Bose exchange symmetry, and the additional Ward 
identity following from general coordinate invariance are sufficient to yield the complete {TJJ) amplitude, 
(4.27), given by eqs. (4.39)-(4.41) with (4.36), (4.33) and Tables IV and|V] for any value of mass or the 
kinematic invariants, without any need of regularization of ultraviolet divergent integrals. 



In Sec. 5 we computed the trace and found the finite anomaly (5.17), equivalent to (4.6). The coefficient 



of the second possible tensor in the trace defined in (4.20) is non- anomalous, but both are needed to 
determine the scaling violation /3 function at finite momentum and finite electron mass. This infrared /? 



function is given by (5.21) in terms of the photon polarization which vanishes when <^ m?, consistent 
with decoupling, and approaches the more commonly considered ultraviolet /3 function only in the opposite 
limit ^ m^. 



In Sec. 6 we gave the spectral representation (6.2) for the imaginary part of the triangle amplitude. 
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cut as in Fig. [2| for = —s < timelike. We showed that the imaginary part of the amphtude is 
non- anomalous, with the anomaly in the real part arising from a cancellation between factors oi k'^ + s in 



both the numerator and denominator of (6.2). For the particular linear combination of spectral functions 



appearing in the anomalous trace, prp{s) defined by (6.4a), we derived on the one hand the finite sum rule 



(6.5), and on the other hand the representation (6.6), which shows that Pj, must develop a 6{s) singularity 



when the fermion mass, and photon virtualities p"^ and q'^ vanish. We also exhibited this 5{s) explicitly in 



this limit, (6.7). Corresponding to the 6{s) in the spectral function p^,, the corresponding full amplitude 
(6.9) has a pole at k'^ = 0, indicating the presence of a massless scalar propagating state in the matrix 



element of the stress tensor to physical photons (6.18). As in the chiral case the massless anomaly pole 



can be understood as arising from a correlated two-particle coUinear e~^e~ state, which because of the 
kinematics is essentially 1 + 1 dimensional, c.f. Fig. [7} 

In Sec. 7 we showed that the trace part of the {TJJ) triangle amplitude is identical with that predicted 



by the covariant effective action (7.1), obtained in earlier work by integrating the anomaly. In particular. 



the variation of the simplified effective action (7.7) in terms of the two scalar auxiliary fields ■0' and (p 



yields the amplitude (7.11) which coincides with the first term of (4.11) which gives its full trace in the 



massless limit. In the effective action the massless scalar two particle state of the triangle amplitude is 
replaced by scalar fields, whose propagator gives rise to the anomaly pole at k'^ = 0. 



Finally in Sec. 8 we considered the tree level gravitational scattering amplitude (8.1), Fig. 9 with one 



vertex replaced by the triangle amplitude {TJJ), and showed that the massless scalar pole in the latter 
survives in the physical scattering amplitude in the conformal limit of massless electrons. In the effective 



theory (8.9) it is described as a propagating massless scalar interaction. Fig. 11 between the trace parts of 



the energy-momentum sources. Abstracting from the axial and trace anomaly QED examples presented 
here in detail, we conclude that the trace anomaly and anomaly pole of conformal fields imply the existence 
of new long range scalar interactions with a gravitational coupling strength to ordinary matter. 
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APPENDIX A: EXTRACTION OF FINITE PARTS OF (TJJ) 



For the amplitude (4.11), in order to extract the finite terms for which each of the indices {^va[3) is 



associated with an external momentum p or q, we may drop the g^^ terms in the vertex y^^, and consider 
only 



dH tr {[Yi'il +p-qY + 7^{2l +p- q^] (-/ - 1^ + m) 7° (-/ + m) 7^ (-/ + ^ + m) } 



(Al) 



4 J (2^)4 [(/+p)2+^2][(/_g)2 + ^2][/2+^2] 

where the continuation to Euclidean / has already been performed. Introducing the Feynman parameteri- 



zation, 



(A2) 



with I' = l+px — qy and D given by (2.20), we shift the integration variable in (lAll) from / to Dropping 



the terms involving either powers of /' or m in the numerator, since these cannot give rise to terms which 
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are homogeneous of degree 4 in p and q in ^Y'^'^'^^, and evaluating the finite Euchdean integral, 

dH' 1 1 1 



(27r)4 (/'2 + /))3 327r2Z) 



(A3) 



we obtain from ( Al |, 

^2 



327r2 



1 



1 f'l — X 

dx dy^ - 2x) - q{l - 2y)] - x) + (72/] a [pa; - qy]p\px + (?(!- 2/)]<xtr{7^)7S"7''7'^7"} • 

JO ^ 

(A4) 



Of the 15 terms in the 7-matrix trace, only the 6 terms. 



up f aX Ba n'^'^ nf^^\ 



■g- - g'^^gn - g^'^g'^'^g"' + 5"V") + • • • (A5) 



need to be retained, since the other 9 contract at least two of the free indices (/Uz/a/3) and do not give rise 



to terms in which all four indices are carried by p or q. Thus we retain from (A4) only the terms, 

o2 



e 
8^ 







dx I 
Jo 






3(1 -2x) -q 


(l-2y)]('^[p(l- 


x) + qyY^\px - 


- qyYipx + q{l - y)] 




-2x) - 




-2y)]("[p(l 


- x) + qyY\px - 


- qyf[px + q{i 


-y)r 


-[p{l 


- 2x) - 


q{l 


-2y)]^^[p{l 


- x) + qyX^lpx - 


qyY\px + q{l 


-y)f 




-2x) - 


9(1 


-2y)]("[p(l 


- x) + qy\^\px - 


qyY\px + q{l 


-y)r 


-[p{l 


-2x) - 


q{l 


-2y)](^[p(l 


- x) + qyX^lpx - 


qyY\px + q{l - 


-y)P 


-[p{l 


-2x) - 


-9(1 


-2y)]('^[p(l 


— x) + qy\^[px — 


qyY[px + q{l 


-yT^} ■ 



(A6) 

In this form it is straightforward to collect the terms which multiply each of the 12 tensors of degree 4 
which are listed in Table [v| of the text. For example, the coefficient of the tensor pi^p^p'^pP from (A6) is 

£^j\xj^' ^dy^ {-4(l-2x)(l-x)x2} . 



(A7) 



When we add the Bose symmetric contribution to (A6) with p replaced by q and a replaced by /?, the 



coefficient of q^q'^q°^ql^ will give an equal contribution to the coefficient Ci, after also interchanging the 
parameter integration variables x and y. Thus, 

ci{x,y) 



Ci{k^;p^q^' 



. o I dx 
47r2 Jo JO 



dy 



D 



(A8) 



with 



ci(x,y) = -4x^(1 -x)(l - 2x) , (A9) 
and we have verified (4.33) for the first entry of Table [V| of the text. The remainder of the Table [V| may 



be derived from (A6) in the same way. 
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For the amplitude A"^ the calculation is similar. Beginning with (5.1), we have 
A°^(p, q) = e^m 



dH tr {(-/ - |5 + m) 7" (-/ + m) -^<^ + )^ + m)} 



+ (p^g,a^/3). (AlO) 



(27r)4 [(/ + p)2 + m2] [(/ - qf + \l? + 

Since the trace of an odd number of 7 matrices vanishes, only those terms with at least one additional 
factor of m in the numerator survive. Since we wish to extract only those finite terms homogeneous of 
degree 2 in the external momenta, namely q'^pp or determining the other terms by vector 



current conservation (5.2 ), we focus only on those terms with exactly one additional factor of m. Using the 



Feynman parameterization (A2), shifting integration variables from Z to /' = Z + — qy^ and evaluating 



the momentum integral (A3) as before, we find from the first term of (AlO), 
dx J dy — |(px - qy)x[px + q{l - y)]/,tr(7"7^7^7^) 



167r2 







-[p(l -x)+ qy]x[px + q{l - y)]Ml^l'"^'^l'") " ^(1 - x) + qy]x{px - gy)^tr(7S"7V)}(All) 
In the 7 matrix traces we may further discard all terms involving g"^, which leaves the remaining terms. 



2 2 />i 



47r2 



dx 



l-x 



dy {{px- qy)"[px + q{l - y)f + {px - qy f[px + q{l - y)Y 



-[p{l -x)+ qyr[px + q{l - y)f + [p{l - x) + qyf[px + q{l - y)]" 
-[p{l -x) + qyripx - qyf - [p{l - x) + qyf{px - gy)"| . 

Adding the Bose symmetrized term with p <^ q and a ^ /?, we obtain for these finite terms, 

l — X -j^ 



(A12) 



2 2 l'^ 
em 



2vr2 7o 



dx 



I dy — |2p>^2;(l - 2x) + p"/(l - 2x)(l - 2y) - q"p^{l - 4xy) + 2g°/y(l - 2y) j(A13) 



From the definitions of Gi and G2 in (5.3) and u°'^{p, q) and w'^^{p, q) in (4.20), it follows that the coefficient 
of —q'^pf^ is 



Gi{k'';p^q^] 
and the coefficient of (p • q)p°'q^ is 

G2{k';p',q') = - 



2 2 fi 



27r2 7o 



dx 



l-x 



dy 



^1 - ^xy) 
D 



(A14) 



2'K^p-q Jo 



dx 







(l-2x)(l-2y) 
ay ^ 



(A15) 



which are Eqs. (5.4) of the text. 
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APPENDIX B: PROOF OF IDENTITIES 



The coefficients of tlie p°^p^ and q'^q^ terms in (A13) apparently do not matcli tliose of w°'^{p,q) witli 
tlie identification of G2 given. Tliis mismatcli is only apparent, because of the identities, 



1 , /-i-^' , 2x(l - 2a;) 
dx ay 
Jo 

1 rl-x 

dx dy 
Jo 



D 

2y(l - 2y) 



dxTdy^'-'^^^'-'y^ 

p-q Jo Jo D 

2 n ri-x q_2x)(l-2y) 



-— — I dx „^ 
D P-Q Jo Jo D 

These identities are most easily proven by considering integrals of the kind. 



dy 



(Bla) 
(Bib) 



I'l I'l-x d\nD f^^'' 1 

dx dy{l-2x)^—= dx dy{l-2x)-- 
Jo Jo dy Jo Jq D 



1 dD 
dy ' 



(B2) 



which on the one hand vanishes, because 



Jo 



dx{l - 2x) InD 



^ , , s , f x(l — x)P + 

dx(l-2x)ln' ^ ' 



0, 



(B3) 



y=l—x 

= I dx[l -2x) [n< 
y=0 Jq [ x{l - x)p^ + 

(1 — 2x) is odd upon reflection about the midpoint of the integral. 



due to the fact that (1 — 2x) — 
a; — > 1 — X, whereas x(l — x) ^ x{l — x) is even; while on the other hand, (B2 ) is equal to 

-I 



dx / dy {1 — 2x) 
Jo 



(1 - 2y)q^ +2xp-q 
D 



(B4) 



Setting this expression to zero and rearranging gives (Bla). The second identity (Bib) is proven in a similar 



manner by exchanging x and y. When these two identities are substituted into (A13), the p'^p^ ,p°'q^ and 



q"q^ terms become proportional to the corresponding terms in w'^^ defined in (4.20), and the coefficient 



function (A15) is obtained in every case. We remark also that despite appearances, G2 has no pole at 



p - q = 0, since the integral in ( A15 ) multiplying it vanishes if we set p - q = in the denominator D, by an 



argument similar to that leading to ( B3 ) . 



The identities (4.32) are proven in a similar way. For example, from Table [V| and (4.33), 



(p ■ q)Ci + q^Ci 



27r2 



dx 







„ , dx xil — x){\ — 2x) 
2vr2 Jo JO 



^-"^ , x(l -x)(l -2x) r „ ^ on 

dy ^ ^ [2xp -q+il- 2y)q^] 

, dlnD 
dy 



e 

'2^ 



dy 

x{l — x)k'^ + m? 



2 fi 

dx xil — x){l — 2x) In > 
[ x(l — x)p^ + 



0, 



(B5) 



for the same reason (B3) vanishes. This proves (4.32a), with (4.32b) proven in exactly the same manner 
after interchanging p^ and g^, and x and y. 
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For the first of identities (4.42), we employ a similar method. In the Feynman parameterized integral 



representation, using (|433j) with Table |V) we find that the linear combination, 

p2 [Ci - 2C2 + C3- 2C7 + 2C78) -2p-q{C7- C9) - (208 - 2C9 + Cio - 2Cu + C; 



12 



is proportional to 



[ dx [ 

Jo Jo 



(1 — X — y) 

dy — {x(x + 3y — 1)(1 — 2x)p'^ — 4xy{x — y)p ■ q — y{3x + y — 1)(1 — 2y)g^} 



D 



^ f^-^ , A , BlnD , ^dlnD 

dx dy [1 — X — y) < x{x + 3y — 1) — y(3x + y — 1) 

Jo I 9x 

l-x 



dy 



f dx [ 

Jo Jo 



dy In D 



d_ 
dx - 



x(l — x — y){l — X — 3y) 



d_ 

dy L' 



y{l-x-y){l-3x-y) ^ =0, (B6) 



which vanishes identically. 



Lastly, the linear combination of terms in (5.11) 



p\Ci - 2C2 + C3 + 2C4 - 4C8) + 2p ■ q{Cs - 2C5 -C7-Cg + do) + q\2Ce - 4Cs + Cw - 2Cu + Cu) , 



after substituting for the Cj from (4.33) with Table [V| is proportional to 



l-x 



1 



dx dy — < X 
Jo ^ 



x(l - x) + y(3y - 2) (1 - 2x) / - 2xy x(l - 2x) + y{l - 2y) 



y(l-y)+x(3x-2) {I - 2y) q 



p-q 



1 rl-y 

dy dx x 
Jo 



x(l-x) + y(3y-2) 



dlnD , /"i-^ , 
V dx dy y 

Jo Jo 



y{l - y) + x{3x - 2) 

[ dyy{l-y){l-2y)ln[y{l-y)k'^ +171^] - [ dx x(l - x)(l - 2x) In [x(l - x)fc^ + 
Jo Jo 



(B7) 
dlnD 



dy 

,21 



mj — / axx(i — x)(i — /x)m|x(i — x)/>; +m 
/o " Jo 

+ dy dx{x-y){3y + 3x-2)lnD- dx dy (x - y){3y + 3x - 2)ln D = , 
Jo Jo Jo Jo 

which also vanishes because the last two terms cancel, while the first and second terms are each separately 

zero by their odd parity under reflection through the midpoint of the remaining integral. The other 

identities can be checked by similar methods, and with the help of algebraic manipulation software such 

as Mathematica. 



